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INTRODUCTION 


The Demonstration Problem Manual is one of four manuals that constitute the documentation for 
NASTRAN. The other three are the Theoretical Manual, the User's Manual and the Programmer's Manual 

The Theoretical Manual contains discussions of the underlying theory relative to the eingineer 
ing equations utilized and mathematical operations. There is some discussion relative to data 
processing techniques and software organization. 

The User's Manual is an instructional and encyclopedic reference that describes finite element 
modeling techniques and shows the requirements for data preparation to obtain solutions in several 
engineering disciplines. 

The Programmer's Manual contains descriptions of the Functional Modules, subroutines and oper- 
ating systems from a software point of view. It also contains detailed derivations of the mathe- 
matical equations employed by the program. 

The Demonstration Problem Manual illustrates the types of problems that can be solved with 
NASTRAN and shows that the results obtained are valid. Generally, this manual discusses the nature 
of the problem, the underlying theory, the specific geometric and physical input quantities, and 
the comparison of theoretical and NASTRAN results. At least one problem for each of the rigid 
formats and nearly all of the elements is provided. The features of NASTRAN demonstrated by speci- 
fic problems are listed in tables which follow. 

The data decks necessary to execute these problems are contained on a tape compatible with 
each of the NASTRAN computers (See the Programmer's Manual, Section 5.4 for system descriptions.) 
The driver decks include the Executive, Case and Substructure Control decks plus changes to the 
Bulk Data (where restarts are involved). The Bulk Data decks are contained on a NASTRAN generated 
UMF (User Master File). To obtain the decks, it is necessary to print the contents of the tape 
8 y using the provided driver decks, NASTRAN may be executed as a UMF 30 b to obtain the results 
for a particular demonstration problem (or the user's desired variation thereof). 

Each demonstration problem is assigned a problem number to key it to the Rigid Format. In 
turn, the UMF problem identification ( pi d) is an adaptation of the problem number The UMF tape 
identification (tid) is the year in which the set of demonstration problems was generated so this 
would tend to change from level to level of NASTRAN caused by the inclusion of new capabilities. 
Furthermore, it may not always be possible to execute the new UMF on a previous level due to 
changes in data handling techniques. 
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A UMF problem identification number is made up of four elements: The Rigid Format number, 

a problem number, the version number, and a trailing dummy zero. Thus the general UMF number 
> xxyyzO. The Rigid Format number is one or two digits; the problem number is always two digits; 
,he version number is always one digit; and the 0 always trails to allow the insertion of additional 
problems. A UMF pi d of 10210 means the problem runs on Rigid Format 1, it is the second demonstra- 
tion problem on that Rigid Format, and it is version 1 of that problem. Another example, 110110, 
is a problem for Rigid Format 11, problem 1, version 1, and the trailing zero is a dummy. 

A table of pi d numbers for each demonstration problem is included in this manual. Restart 
problem driver decks do not contain a UMF card because the data is already stored on a checkpoint 
tape which must have been created by the user. 
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Windmill Panel Sections for Multi-stage Substructuring (Run 6, Phase 3) 

Windmill Panel Sections for Multi-stage Substructuring (Run 7, Phase 2) 

Vibration of a 10x20 Plate 

Vibration of a 20x40 Plate 

Vibration of a 10x20 Plate (via INPUT Module) 

Vibration of a 20x40 Plate (via INPUT Module) 

Vibration of a Compressible Gas in a Rigid Spherical Tank 

Vibration of a Liquid in a Half Filled Rigid Sphere 

Acoustic Cavity Analysis 

Nonlinear Heat Transfer in an Infinite Slab 
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DEMONSTRATED FEATURES OF NASTRAN 


A. PHYSICAL PROBLEMS 
Structures 


1 . 

Line 

2. 

Plate or Shell 

3. 

Solids 

4. 

Rotational Symmetry 

Fluid Dynamics 

5. 

Flexible Boundary 

6. 

Rigid Boundary 

7. 

Sloshing 

8. 

Acoustic 

9. 

Aeroelastic 

Heat 

Transfer 

10 . 

Conduction 

n. 

Convection 

12. 

Radiation 

SOLUTION METHODS 

Steady State 

1 . 

Linear Statics 

2. 

Inertia Relief 

3. 

Nonlinear Geometry 

4. 

Material Plasticity 

5. 

Fully Stressed Design 

6 . 

Linear Heat Transfer 

7. 

Nonlinear Heat Transfer 
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8. 

Real Modes 

9. 

Complex Modes 

10 

Inverse Power 

11. 

Determinant 

12. 

Givens 

13. 

Upper Hessenberg 

Dynamic Response 

14. 

Direct Formulation 

15. 

Modal Formulation 

16. 

Transient Response 

17. 

Frequency Response 

18. 

Random Analysis 

19. 

Flutter Analysis 

ELEMENT TYPES 

1 . 

Bar, Rod, Tube or Conrod 

2. 

Shear or Twist Panel 

3. 

Plate or Membrane 

4. 

Scalar Springs, Mass and Dampers 

5. 

Concentrated Mass 

6. 

Viscous Dampers 

7. 

Plot (PL0TEL) 

8. 

General (GENEL) 

9. 

Conical Shell 

10. 

Toroidal Shell 

11 . 

Axi symmetric Solids 

12. 

Linear Sol ids 

13. 

Isoparametric Solids 

14. 

Solid Heat Conductors 

15. 

Heat Transfer Boundary Elements 
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demonstrated features of nastran 


16. Fluid Elements 

17. Acoustic Elements 

18. Aerodynamic 

D. CONSTRAINTS 

1. Single-Point Constraints 

2. Multipoint Constraints 

3. Omitted Coordinates 

4. Free-Body Supports 

5. Fluid Free Surface 

6. Symmetry Used on Boundary 

7. "Grounded" Stiffness Terms 

E. GEOMETRY AND PROPERTY DEFINITIONS 

1. Property ID Default 

2. Local Coordinate System 

3. Resequenced Grid Points 

4. Thermal Dependent Materials 

5. Nonlinear Materials 

6. Anisotropic Materials 

7. Offset BAR Connections 

8. Structural Mass 

9. Nonstructural Mass 

10. Structural Element Damping 

11. Compressibi 1 ity of Fluid 

12. Fluid Gravity Effects 

13. Multiple Fluid Harmonics 


6 (3/1/76) 



DEMONSTRATED FEATURES OF NASTRAN 


F. SPECIAL MATRIX OPTIONS 

1. General Element (GENEL) 

2 Direct Input Matrices 

3. Transfer Functions 

4. Extra Points 

5. Direct Damping Matrix Input 

6. Modal Damping 

7. Substructuring 
8 Cyclic Symmetry 

9. Uniform Structure Damping 

G. LOADING OPTIONS 
Static 

1. Concentrated Loads 

2. Pressure Loads 

3. Gravity Loads 

4. Thermal Loads 

5. Harmonic Loads 

6. Centrifugal Field Loads 

7. Enforced Element Deformation 

8. Enforced Displacement 

Dynamic Excitation 

9. Tabular Loads vs. Frequency or Time 

10. Direct Time Function Loads 

11. Loading Phase Angles 

12. Loading Time Lags 

13. Load Combinations (DL0AD) 

14. Transient Initial Conditions 

15. Random Analysis Power Spectral Density Functions 
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DEMONSTRATED FEATURES OF NASTRAN 


Heat Transfer 

16. Volume Heating 

17. Area Heating 

18. Radiation Heating 

19. Enforced Boundary Temperature 

H. EXECUTION OPTIONS 

Multiple Solution Techniques 

1 . Loads 

2. Boundary Constraints 

3. Cyclic Cymmetry 

4. Direct Input Matrices 

5. Aerodynamic Coefficients 

Operational Techniques 

6. Checkpoint 

7 Restart with Modified Case Control 

8. Restart with Rigid Format Change 

9. Restart with Modified Bulk Data 

10. Altered Rigid Format Using DMAP Statements 

11 . Multi-stage Substructuring 

I. OUTPUT OPTIONS 
Print an d/o r Punched 

1. Point Output Selections 

2. Element Output Selections 

3. Subcase Level Request Changes 

4. Sorted by Frequency or Time (S0RT2) 

5. Magnitude and Phase of Complex Numbers 

6. Mode Acceleration Data Recovery 
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DEMONSTRATED FEATURES OF NASTRAN 


7. Solution Set Output Requests 

8. Frequency Set Selections 

9. Punched Output Selections 

10. Weight and Balance 

11. Grid Point Force Balance 

12. Element Strain Energy 


Plot 


13. 

Structures Plot of 

Undeformed Structure 

14. 

Structures Plot of 

Deformed Structure 

15. 

Curve Plotting vs. 

Frequency 

16. 

Curve Plotting vs. 

Time 

17. 

Curve Plotting vs. 

Subcase 
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Note: Two-digit problem numbers refer to all problems in the series: three-digit numbers refer 

to a specific version. 
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modal formulation 


Note: Two-digit problem numbers refer to all problems in the series; three-digit numbers 

refer to a specific version. 
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Note: Two-digit problem numbers refer to all problems in the series; three-digit numbers refer to a specific version. 
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Note: Two-digit problem numbers refer to all problems in the series; three-digit numbers refer to a specific version. 
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Note: Two-digit problem numbers refer to all problems in the series; three-digit numbers refer to a specific version. 
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Note: Two-digit problem numbers refer to all problems in the series; three- 

digit numbers refer to a specific version. 
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Note: Two-digit problem numbers refer to all problems in the series; three-digit numbers refer to a specific version. 
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numbers refer to a specific version. 
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Note: Two-digit problem numbers refer to all problems in the series; three-digit numbers refer to a specific version. 
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RIGID FORMAT No. 1, Static Analysis 

Delta Wing with Biconvex Cross Section (1-1-1) 

Delta Wing with Biconvex Cross Section Using QDMEM1 and QDMEM2 Elements (1-1-2) 

Delta Wing with Biconvex Cross Section Using QDMEM1 Elements (1-1-3) 

Delta Wing with Biconvex Cross Section Using QDMEM2 Elements (1-1-4) 

This problem illustrates the use of various NASTRAN elements in the solution of an actual 
structural problem. Figure 1 shows the delta wing to be modeled and Fiaures 2 and 3 shows the 
finite element model. The delta wing model is composed of membrane, shear panel and rod elements. 
Due to the existence of symmetry or antisymmetry in the structure and loadinq conditions, only 
one-quarter of the wing needs to be modeled. The midplane of the wing (the plane dividing the wing 
into upper and lower halves) is a plane of symmetry as is the center plane (the plane that divides 
the wing into left and right halves). The loading conditions are anti symmetrical with respect to 
the midplane of the wing and symmetric with respect to the center plane. 

The surface skin of the wing is modeled with membrane elements while the ribs and spars are 
modeled with a combination of shear panels and rods. The shear load carrying capability of ribs 
and spars is represented by shear panels. The bending stiffness of the ribs and spars is modeled 
with rod elements placed in the plane of the skin surface. 

Since a quarter model is used, the loading conditions require that an antisymmetri c boundary 
be provided on the midplane and a symmetric boundary must be provided on the center plane. These 
boundary conditions are provided by constraining all grid points on the midplane in the x and y 
directions and all grid points on the center plane in the x direction. Supports for the structure 
are provided by constraining grid points 13, 33, 53, 73 and 93 in the z direction. Since no rota- 
tional rigidity is provided by the elements used in the model, all rotational degrees of freedom 
have been removed by the use of the GRDSET card. 

Figure 4 shows the two loading conditions analyzed. The modified restart capability is used 
to perform the analysis associated with the second loading condition. The ability of NASTRAN to 
change rigid formats on a restart is demonstrated by the third case. The natural modes of the 
structure are extracted using the inverse power method. Since the symmetric boundary conditions 
are used, only the modes with symmetric motion about the center line will be extracted. If the 
unsymmetric modes were required, a separate run with the aporopriate boundary conditions could be 
submitted. 
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A comparison of the displacements due to the loads calculated by NASTRAN and the experimentally 
measured displacements from Reference 2 is shown in Tables 1 and 2. Also included in these tables 
are the displacements calculated on a passive analog computer (Reference 1). A comparison of the 
natural frequencies and modal displacements is shown in Tables 3 and 4. 

Another variation of this problem can be obtained by replacing the quadri lateral membrane 
elements (CQDMEM) with the newer CQDMEM1 and CQDMEM2 elements. This modification demonstrates the 
ability to reproduce previously derived theoretical results. Table 5 shows the difference in dis- 
placements obtained when elements 1 through 9 are CQDMEM1 elements and the other quadrilaterals 
are CQDMEM2 elements. 
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Table 1 Comparison of NASTRAN and Experimental Deflections - Concentrated Load on Outboard 

Trailing Edge 


GRID 

NUMBER 

Z DISPLACEMENT 

NASTRAN 

EXPERIMENTAL 

ANALOG 

14 

-.082 

-.08 

-.080 

15 

- 221 

-.22 

-.210 

16 

-.424 

-.39 

-.400 

34 

'-.063 

-.07 

-.061 

35 

-.162 

-.16 

- 157 

36 

-.293 

-.28 

-.286 

54 

-.043 

-.05 

-.044 

55 

-.104 

-.12 

-.144 

74 

-.025 

-.03 

l 

-.030 


Table 2 Comparison of NASTRAN and Experimental Deflections - Concentrated Load on Outboard 

Leading Edge 


GRID 

NUMBER 

Z DISPLACEMENT 

NASTRAN 

EXPERIMENTAL 

ANALOG 

14 

-.063 

- 06 

- 060 

15 

-.163 

-.15 

- 157 

16 

- 293 

-.28 

- 286 

34 

- 057 

-.06 

- 057 

35 

-.148 

-.15 

-.150 

36 

-.280 

- 30 

-.290 

54 

-.046 

-.05 

- 048 

55 

- 118 

- 13 

-.127 

74 

- 030 

-.04 

- 035 
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TABLE 3 Comparison of NASTRAN and analog computer analysis eigenvalues. 


Mode No. 

NASTRAN (cps.) 

ANALOG (cps.) 

1 

40.9 

WBBKM 

2 

115.3 


3 

156.2 

' 1 


TABLE 4 Comparison of mode displacements for first mode. 


GRID 

NUMBER 

Z DISPLACEMENT 

NASTRAN 

ANALOG 

14 


.273 

15 


.630 

16 


1.000 

34 

.210 

.239 

35 

.504 

.558 

36 

.854 

.902 

54 

.162 

.192 

55 

.391 

.462 

74 

.112 

.148 
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TABLE 5. Comparison of Z Displacements 


Trailing Edge Load 

Leading Edge Load 

CQDMEM Elements 

CQDMEM1 and 
CQDMEM2 
Elements 

CQDMEM Elements 

CQDMEM1 and 
CQDMEM2 
Elements 

- * 082 

- 082 

- 063 

-.064 

-.221 

- 224 

-.163 

-.167 

-.424 

-.433 

- 293 

-.300 

-.063 

-.064 

- 057 

-.059 

-.162 

- 166 

- 148 

-.154 

-.293 

- 300 

-.280 

-.294 

-.043 

- 044 

- 046 

- 047 

-.104 

-.108 

-.118 

- 123 

-.025 

-.026 

- 030 

-.031 
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Figure 1. Delta wing with biconvex section 
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Delta wing with biconvex section model. 
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RIGID FORMAT No. 1, Static Analysis 
Spherical Shell with Pressure Loading (1-2-1) 


This problem demonstrates the finite element approach to the modeling of a uniform spherical 
shell. A spherical coordinate system is chosen to describe the location and disnlacement deqrees 
of freedom at the grid points. Triangular plate elements are chosen to provide a nearly uniform 
pattern. Two symmetric boundaries are used to analyze the structure with a symmetric pressure 
load. Figure 1 describes the model. 

Two boundary conditions are used on the outside edge to demonstrate the ability of NASTRAN 
to restart with different constraint sets by simply changing the case control request. The 
effective boundary constraints are shown in Figure 2. The membrane support, under a uniform 
pressure load, should result in uniform in-plane compression in two directions. The clamped 
support produces bending moments in addition to in-plane stresses. 

The grid point numbering sequence used minimizes the computer time required to perform the 
triangular decomposition of the constrained stiffness matrix. This numbering sequence results in 
a partially banded matrix with all terms outside the band located in a single column. The grid 
points are arranaed to form five rings; the center point is sequenced last. 

Analytic solutions for the continuum shell were obtained from Reference 4 using the first 20 
terms of the series shown in Equation (j ) of Section 94. Comparisons of the answers obtained 
using NASTRAN and the analytical solution for the membrane boundary condition are shown in Figures 
3 and 4. Also included on these figures are the NASTRAN answers obtained usina a 10-ring model 
Figures 5 thru 7 show a comparison of the NASTRAN answers and the analytical solution for the 
shell with a clamped boundary. 

The slight differences between theoretical and computed answers are due to* 1) The finite 
element model assumes a constant in-plane stress and a linearly varying bending moment for each 
element. In the clamped edge case these quantities have large chanaes , and 2) the irregular! ties 
of the finite element model cause some extra couplinq between the bendinq and membrane action. 
Since the elements are planar the curvature is modeled, in effect, by the dihedral angles between 
elements. Since the elements are different sizes and shapes these dihedral anales vary, which 
results in slight differences in curvature that cause small errors. 
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Figure 1. 5 ring spherical shell model. 
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Load 1 lb/in.2 



(Problem l-2a) (Problem 1-2) 


Figure 2. Spherical shell loading and edge suoport conditions. 
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Figure 5. Comparison of NASTRAN and analytical displacements for 
5 ring spherical shell - clamped boundary. 
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Figure 6. Comparison of NASTRAN and analytical meridian stress for 
5 ring spherical shell - clamped boundary. 
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Figure 7. Comparison of NASTRAN and analytical bending moment for 
5 ring spherical shell - clamped boundary. 


1.2-7 (6/1/72) 



RIGID FORMAT No. 1, Static Analysis 

Free Rectangular (QDMEM) Plate with Thermal Loading (1-3-1) 

Free Rectangular (QDMEM1) Plate with Thermal Loading (1-3-2) 

Free Rectangular (QDMEM2) Plate with Thermal Loading (1-3-3) 

This problem demonstrates the use of thermal loading conditions and temperature dependent 
materials. The model, a rectangular plate shown in Figure 1, is given a temperature gradient 
which causes internal loads and elastic deflections. Since there are two planes of symmetry, 
only one-quarter of the structure needs to be modeled (the shaded portion shown in Figure 1). 

The finite element model for the quarter section is shown in Figure 2. Fiaure 3 shows the 
thermal loading condition. The temperature load is constant in the y direction and symmetric 
about the y-axis. Since membrane elements are used to model the structure, it is necessary to 
remove all rotational degrees of freedom and translational degrees of freedom normal to the plate. 
The symmetric boundary conditions were modeled by constraining the displacements normal to the 
planes of symmetry. The material used has temperature dependent elasticity (as defined in Reference 
5) therefore the INPUT module cannot be used for this application. The CNGRNT bulk data card can 
be used if the congruency is defined in one direction. 

Figures 4 and 5 show a comparison of NASTRAN stresses and the experimentally measured stresses 
reported in Reference 5. 
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Figure 1. Free plate structure 
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Figure 4. 


Comparison of NASTRAN and experimental stresses for free rectangular 
plate with thermal loading - temperature dependent properties. 
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Figure 5. Comparison of NASTRAN and experimental stresses for free rectangular 
plate with thermal loading - temperature dependent properties 
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RIGID FORMAT No. 1, Static Analysis 

Long, Narrow, 5 x 50 Orthotropic Plate (1-4-1) 

Long, Narrow, 5 x 60 Orthotropic Plate (1-4-2; 

Long, Narrow, 5 x 50 Orthotropic Plate (INPUT, 1-4-3) 

Long, Narrow, 5 x 60 Orthotropic Plate (INPUT, 1-4-4) 

This problem demonstrates triangular decomposition spill logic, orthotropic materials, and 
the use of a modified restart to obtain additional output. A sketch of the rectangular plate and 
the applied loading is shown in Figure 1. 

The 5 x 50 finite element quarter-section model shown in Figure 2 is constructed with quadri- 
lateral bending elements. In order to demonstrate the triangular decomposition spill logic (i.e. 
the necessary computation space is larger than available core storage), the model is internally 
resequenced for a wide bard as shown in Figure 3. Although the 5 x 50 model is sufficient to 
create spill on the IBM 7094 DCS, the number of elements in the longitudinal direction must be 

increased to create spill on machines having larger random-access memories. 

The analytical solution for the infinitely long continuum plate is given in Section 37 of 
Reference 4. A comparison with the NASTRAN solutions are given in Tables 1 and 2. 

A modified restart was used to obtain additional output. 
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Table 1. Comparison of NASTRAN and Analytical Displacements for Long, Narrow, Grthotropic Plate 


GRID 

Z DISPLACEMENT x 10 4 | 

THEORY 

NASTRAN 

1 

3.048 

3.037 

2 

2.899 

2.889 

3 

2.466 

2.457 

4 

1 .792 

1.785 

5 

0.942 

0.939 

7 

2.949 

2.940 

13 

2.723 

2.714 

19 

2.446 

2.435 

25 

2.157 

2.145 

31 

1 .880 

1.866 

37 

1.625 

1 .611 

43 

1.397 

1.383 


Table 2. Comparison of NASTRAN and Analytical Displacements for Long, Narrow, Orthotropic Plate 


EL. 

ID. 

STRESS X 

STRESS Y 

SHEAR 

STRESS 

THEORY 

NASTRAN 

THEORY 

NASTRAN 

THEORY 

NASTRAN 

1 

19.05 

18.90 

20.35 

20.40 


-0.39 

2 

17,19 

17.05 

18.36 

18.40 

ns 

-1.13 

3 

13.64 

13.53 

14.57 

14.60 

-1.74 

-1.76 

4 

8.76 

8.69 

9.35 

9.38 

-2.19 

-2.22 

5 

3.02 

2.99 

3.22 

3.23 

-2.43 

-2.46 

7 

15.86 

15.76 

12.91 

12.90 

-0.84 

-0.88 

13 

13.27 

13.20 

8.28 

8.23 

-1.03 

-1.06 

19 

11.14 

11.08 

5.38 

5.33 

-1.07 

-1.09 

25 

9.37 

9.33 

3.55 

3.51 

-1.02 

-1.04 

31 

7.90 

7.86 

2.38 

2.36 

-0.94 

-0.95 

37 

6.67 

6.63 

1.64 

1.63 

-0.84 

-0.85 
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Figure 1. Simply-supported long narrow orthotropic plate 
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RIGID FORMAT No. 1, Static Analysis 
Nonsymmetric Bending of a Cylinder of Revolution (1-5-1) 


A. Description 

This problem illustrates the application of the conical shell element and its related special 
data. This element uses the Fourier comporents of displacement around an axisymmetric structure 
as the solution coordinates. The geometry of the structure is defined by rings instead of grid 
points. Its constraints must be defined by the particular Fourier harmonics, and the loads must 
be defined either with special data or in a harmonic form. This element may not be used in conjunc- 
tion with any of the other structural elements. 

The structure to be solved is described in Reference 6 and illustrated in Figure 1. It con- 
sists of a short, wide cylinder with a moderate thickness ratio. The applied loads and the out- 
put stresses are pure uncoupled harmonics. The basic purpose of this problem is to check the 
harmonic deflections, element stresses, and forces. Figures 2 and 3 compare the NASTRAN results 
with the results given in Reference 6. 

B. Input 

The Fourier coefficients of the applied moment per length are: 

m n = cos(ne) 

I 

The applied input loads are defined as: 

2tt 

M n = / m cos(ne) R de 
o 

The values of applied moment on the M0MAX cards are: 


M n0> = 

ttR 

, n > 0 

II 

-e- 

o 

IE 

2irR 

, n = 0 


The applied moments for each harmonic are shown in Figure 1. The bendina moments in the 
elements are defined as: 
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Moment about u 


M = Moment about u 
u z 

Positive bending moments indicate compression on the outer side. 

1. Parameters: 

R = 50 Radius 

s = 50 Height 

t = 1.0 Thickness 

E = 91.0 Modulus of Elasticity 

v = 0.3 Poisson's Ratio 

2. Loads: 

M n(100) = 157.0796 Force -Length 

M n(50) ~ Force-Length 

3. Single Point Constraints: 


Ring ID 

Harmonic 

Coordinates 




50 

all 

vv u z 

Radial, tangential 

and axial 

translations 

100 

all 

W u z 

Radial, tangential 

and axial 

translations 

all 

all 

9 r 

Rotation normal to 

surface 



The AXISYM = C0SINE statement in case control defines the motions to be symmetric with 
respect to the x-z plane. 

C. Answers 

Theoretical and NASTRAN results for element bending moments and radial deflections for 4 of 
the 20 harmonics used are given in Figure 2 and 3. Notice that for hiqher harmonics the effect of 
the load is limited to the edges. A smaller element size at the edges and a relatively large 
size in the center would have given the same accuracy with fewer degrees of freedom. 
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Figure 1. Cylinder under harmonic loads. 
1 . 5-3 



Figure 2. 


Element bending moments and radial deflections 
along length of cylinder. 
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Figure 3. 


Element bending moments and radial deflections 
along length of cylinder. 
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RIGID FORMAT No. 1, Static Analysis 
Solid Disk with Radially Varying Thermal Load (1-6-1) 


A solid free circular disk in a plane is subjected to a radially varying thermal load of the 

form 

T = 100(1 - {£) 

where 

r = the radius at any point in the disk 
b = the outside radius = 0 10 inches 


The structure is shown in Figure 1 along with its associated material properties and pertin- 
ent dimensions The finite element idealization employed for this structure is shown in Figure 2. 
The thermal loading on the solid disk is established via an internally generated thermal load 
vector derived from data specified grid point temperature values. 

Figure 3 displays the radial displacement utilizing the idealization shown in Figure 2. 

Figure 4 presents radial and ci rcumferential stress values which result from the thermal loading. 
Reference 14 provides an analytical solution to this problem which is based on the theory of 
elasticity. Note that the solid lines represent the analytical solution while the circles and 
squares represent the solution obtained utilizing the finite element solution 


\ 
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Figure 1. Solid circular disk 


Figure 2. Finite element idealization and 
temperature distribution. 
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Radial Displacement - in. 










RIGID FORMAT No. 1, Static Analysis 

Shallow Spherical Shell Subjected to External Pressure Loading (1-7-1) 

A shallow spherical shell with a built-in edge is subjected to an external pressure loading of 
1 psi . The shell is shown in Figure 1 along with its pertinent dimensions and associated material 
properties. The finite element idealization for the shell is displayed in Figure 2. 

Due to symmetry only one half of the shell was analyzed. The primary purpose of this analysis 
was to demonstrate the applicability of the shell cap generalization of the toroidal ring to this 
class of problem. 

The meridional bending moment is taken to characterize the behavior predicted for this struc- 
ture. The exact solution from Reference 4 is compared to the 13 element finite element solution 
in Figure 3. The reference solution is designated by the solid line while the finite element sol- 
ution is designated by the circles. Figure 4 displays the radial displacement obtained utilizing 
this idealization and compares it to that obtained in the Reference Solution. 
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Figure 4. Radial displacement, shallow 
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RIGID FORMAT No. 1, Static Analysis 
Bending of a Beam Fabricated from HEXA1 Solid Elements (1-8-1) 


A. Description 

The properties of solid bodies may be modeled with the NASTRAN tetrahedra, wedge, or hexa- 
hedron finite elements. This problem demonstrates the analysis of a solid fabricated from the 
six-sided HEXA1 solid elements. The problem consists of a rectangular parallelepiped subdivided 
into forty cubic subelements as shown in Figure 1. 

The loads were chosen to approximate the stress distribution due to a moment on one end of 
a beam; the other end is constrained to resist the moment. Two planes of symmetry were used to 
simulate an actual problem having twice the width and twice the height. 

B. Input 

1. Parameters: 

l = 20.0 (length) 
w = 4.0 (width of full section) 
h = 16.0 (height of full section) 

E = 3.0 x 10 6 (modulus of elasticity) 
v = 0.2 (Poisson's ratio) 

2. Boundary Constraints: 

on y = 0 plane, u x = u z = 0 (antisymmetry) 

on z = 0 plane, u 2 = 0 (symmetry) 

on x = 0 plane, u x = 0 (symmetry) 

3. Loads: 

Total Moment: M y = 2.048 x 10 3 

This moment will produce bending about the z axis. It is modeled by a set of axial 
loads at x = i which, in turn, represent an axial stress distribution: 

a xx = 1 - 5 y 

C. Analysis and Results 

A prismatic beam with an axial stress which varies linearly over the cross section has an 
exact solution. In the demonstration problem, the theoretical stress distribution is: 
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a xx " ' I y 


a yz a zz T xy T xz T yz ® ' 


1 3 

where I = wh • 


The displacements are: 


M 

n- xy 


u y = 2n"(x 2 -vy 2 -vz 2 ) 


M 

u z = et yz 


Tables 1 and 2 are comparisons of displacements and stresses for the theoretical case and 
the NASTRAN model . 
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DISPLACEMENT x 10" 4 

POINT/DIRECTION 

THEORY 

NASTRAN 

21 /y 

.0400 

.0417 

41 /y 

.1600 

.1607 

61 /y 

.360 

.366 

81 /y 

.640 

.651 

101/y 

1.000 

1.016 

109/x 

0.800 

0.844 

110/z 

.016 

0.007 


Table 1. Comparisons of Displacement 


ELEMENT 

°xx 

THEORY 

a yy 

T 

xy 

' a xx 

NASTRAN 

a yy ' 

H 

X 

.*< 

1 

-1.5 

0 

0 

-1.56 

.02 

. -.01 

2 

-4.5 

0 

0 

-4.53 

.036 

-.05 

3 

-7.5 

0 

0 

-7.39 

.06 

-.06 

4 

-10.5 

0 

0 

-9.95 

'-.11 

.12 


NOTE: NASTRAN stresses are average; theoretical stresses 

are calculated at the center of the element. 


Table 2. Comparisons of Stress 
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RIGID FORMAT No. 1, Static Analysis 
Thermal and Applied Loads on HEXA2 Solid Elements (1-9-1) 


A. Descri ption 

This problem demonstrates the use of the HEXA2 solid hexahedron elements. Forty rectangular 
elements are used to model a 2x2x10 beam. The dimensions and boundary conditions are shown in 
Figure 1. Two loading conditions are applied: axial stress and thermal expansion. Symmetry 

boundary conditions are used. 


B. Input 

1. Parameters: 

l = 20 (length) 

.w = 4.0 (width) 

h = 4.0 (height) 

E = 3.0 x 10 6 (modulus of elasticity) 

v = 0.2 (Poisson's ratio) 

a = .001 (thermal expansion coefficient) 

T 0 = 10° (reference temperature) 

2. Boundary Constraints: 
u x = 0 at x = 0 

Uy = 0 at y = 0 
u z = 0 at z = 0 

3. Loads: 

Subcase 1 , 

F x = 24 x 10^ (total axial force) 
Subcase 2, 

T = 60° (uniform temperature field) ' 
T o = 10° (reference temperature) 


1.9-1 (3/1/76) 



C. Analysis and Results 


1 . Subcase 1 

The distributed axial load is equivalent to a stress field of: 


°xx = ] - 5xl ° . 


a = a = t = t_ = t = 0 

yy ZZ xy xz yz 


The displacement field is: 


"x ■ f' 


0.5 x 10 _3 x 


-vo 




-0.1 x 10’ 3 y 


-vcr 


xx 


z = 


-0.1 x 10" 3 z 


2. Subcase 2 

The uniform expansion due to temperature will not cause any stress. The strains, however, 
are uniform and equal. The displacements are, therefore: 

u x = a(T-T Q )x = .05x 


U y = «(T-T 0 )y = -05y 


u z = a(T-T Q )z = .05z 


where T is the uniform temperature and T q is the reference temperature. 


3. Results 

The results of both subcases are exact to the single precision limits of the particular 
computer used. 
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Figure 1 Model of solid using HEXA2 elements 
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RIGID FORMAT No. 1, Static Analysis 
Thermal Bending of a Beam (1-10-1) 


A. Description 

This problem demonstrates the solution of a beam subjected to a thermal gradient over the 
cross-section. Two end conditions are solved, clamped-free and cl amped-pinned end conditions. 

An equivalent linear gradient in the normal direction was used for the input data. However, 
the actual temperatures at points on the cross-section were input on the TEMPRB card in order to 
produce correct stresses. The beam was subdivided into 14 variable lengths for maximum efficiency. 

B. Input 

Figure 1 describes the beam and the thermal field to be analyzed and Figure 2 shows the finite 
element model . 


C. Theory 

For subcase 1, the effective temperature gradient, T' , (see NASTRAN Theoretical Manual) is: 


where 


T'(x) 


} | | T(x,y,z)y dy dz 

z y 


I 



dy dz 


z y 


0 ) 

( 2 ) 


Using the given temperature distribution the effective gradient is: 


T' = T c X 3 , (3) 

where T c is calculated to be 0.170054°/in 4 by substituting the temperature distribution into 
Equation 1 and evaluating the expression: 


T 


c 


1 

I 


j Cy 4 dy 
z y 


dz 


(4) 
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Since the bar is not redundantly constrained the curvature at the center line is: 


d 2 v 

dx 2 


= -aT‘ 


-aT x" 
c 


( 5 ) 


The slope is: 


dv 

dx 


d 2 v 

dx 2 


dx = 


_ a t q 
4 'c X 


( 6 ) 


The deflection is: 


v(x) = 


dv dx = 
dx ax 


- — T x b 
20 'c X 


(7) 


The moment, M, shear, V, and axial stress, a x , are: 


■ EI (^? + “ r ) ' 0 ' 


« ■ f ■ 0 


o x (x,y) = E(£ x ■ aT) = E(ayT ' " aT) = Ect(T c y ■ Cy3) x3 


( 8 ) 


where C = 1 has dimensions of degrees/length . 

For subcase 2, with a simple support at x = 10.0, we calculate the deflection due to subcase 1 
and apply a constraint load to remove the deflection at the end. 

aT 


3EI 


v (L ) 


3EI — 9. j_2 

2Q L 


(9) 


Note: Transverse shear deflection is neglected. 

The deflections and slopes are the sum of the results for the two independent loads as follows 


deflection: 


"lope: 


v(x) = (3Lx 2 - x 3 ) - ^ x 5 = ^ (3L 3 - L 2 x - 2x 3 ) x 2 , 


0.(x) = |r = ur - 3L^x - 10x J )x 


aT. 


3x 


40 


( 10 ) 

OD 
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The net stress is the sum of the stress due to each load: 


° x (x,y) = Ea(T c y - Cy 3 )x 3 - = Ea [( T c y - Cy 3 )x 3 - ± T c L 2 (L - x) y] (12) 

where is the moment due to the constraint load. 

0. Results 

Tables 1 and 2 compare the analytical maximum value of displacement, constraint force, ele- 
ment force, and stress to the maximum deviation of NASTRAN in each category. All results are 
within 2.66%. 
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t 


LULL 


where C = 


y 

I v(x) 


i 


0.2 in. 

pJ i 



1 



n 



0.2 in. -» «- $0.6 in 


1 - - 1 

1 L 1 

e ' 

E = 10.0 x 10 Ib/in (Young's modulus) 0.2 in. 


v = 0.3 (Poisson's ratio) 

a = 13.0 x 10 6 in/1n/°F (Thermal exoansion coefficient) 

Tq = 0.0°F (Reference temDerature) 

The beam is loaded by the temperature distribution: 


1.0 °F/in 6 


T(°F) = Cx 3 y 3 



Figure 1. 

Thermal 

loading of a beam. 
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Figure 2. Finite element model. 
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Table 1. Comparison of NASTRAN and analytical results, clamped-free ends (subcase 1). 


CATEGORY 

MAXIMUM 

ANALYTICAL 

VALUE 

MAXIMUM 

NASTRAN 

DIFFERENCE 

PER CENT 
ERROR 

Displacement 

-1.1054 x 10' 2 

2.9424 x 10~ 4 

2.66 

Constraint Force 

0 

★ 

★ 

Element Force 

0 

* 

-k 

Element Stress 

5.1965 x 10 +3 

0.671 

0.01 


*These results vary with the computer. The very small numbers are essentially zero when 
compared to subcase 2 results. 


Table 2. Comparison of NASTRAN and analytical results, clamped-pinned ends (subcase 2). 


CATEGORY 

MAXIMUM 

ANALYTICAL 

VALUE 

MAXIMUM 

NASTRAN 

DIFFERENCE 

PER CENT 
ERROR 

Displacement 

4.3936 x 10" 3 

8.024 x 10' 6 

0.18 

Constraint Force 

-2.2859 x 10 +2 

6.0841 

2.66 

Element Force 

2.2859 x 10 +2 

6.0846 

2.66 

Element Stress 

5.1965 x 10 +3 

4.4136 x 10 

0.85 
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RIGID FORMAT No. 1, Static Analysis 

Simply-Supported Rectangular Plate with a Thermal Gradient (1-11-1) 

Simply-Supported Rectangular Plate with a Thermal Gradient (INPUT, 1-11-2) 

A. Description 

This problem illustrates the solution of a general thermal load on a plate with the use of an 
equivalent linear thermal gradient. The thermal field is a function of three dimensions, demonstra- 
ted by the TEMPP1 card. The plate is modeled with the general quadrilateral, QUAD1 , elements as 
shown in Figure 1. Two planes of symmetry are used. This problem is repeated via the INPUT module 
to generate the QUAD1 elements. 

B. Input 


E 

= 

3.0 x 10 5 

o 

pounds/inch (Youngs modulus) 

V 

= 

0.3 

(Poisson's ratio) 

P 

= 

1 .0 pound- 

sec. 2 /inch 4 (Mass density) 

a 

= 

0.01 inch/ 

°F/inch (Thermal expansion coefficient) 

t r 

= 

0.0 °F 

(Reference temperature) 

T o 

= 

2.5 °F 

(Temperature difference) 

a 

- 

10.0 inch 

(Width) 

b 

= 

20.0 inch 

(Length) 

t 

= 

0.5 inch 

(Thickness) 

The thermal 

field is: 

' * 




T = T 0 (cos f) (cos f) (ff 

V V 
1 




» 160.0(cos ( C0S IJ) z 3 °F 

s •**, 
] % 

C. Theory 





The plate was solved using a minimum energy solution. The net moments, { } » in the plate 
are equal to the sum of the elastic moments, {M e }, and the thermal moments, } . 

{M n } = {M t } + {M e } , (1) 


where the thermal moment is 
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{M t } = a"P D ( 1 +v ) j 1 1 cos —■ cos ^ 


( 2 ) 


and 


D = 


Et 


12(lV) 


and Tq = 6T Q /5t is the effective thermal gradient. 

The elastic moment is defined by the curvatures, x, with the eauation: 


{M p } = 


*x + vx v 


D < X y + vx x 


(1-v) 


2 X 


xy 


(3) 


Assuming a normal displacement function, W, of 


W = !y T W cos ^x cos m^ 

L L nm a n 


n m 


(4) 


then 


3 2 W 

X x — Z 

X 3x^ 


- 1 1 (?) "P «* T 


n m 


2 

aw 

Xu " 2 


n m 


v r 2,, / m\ nTTX nmy 

* r W nmV?) cos — cos b 


(5) 


X 


= 2 - 


3 2 W 


xy 3x9y 


0 v r 2,, /nm\ nirx „ miry 

2 11 * W n m Uh Sln T- Sln h 


n m i 


The work done by the thermal load is: 


U = 


{ X } T {M t } dA + ij { X } T (M e } 


dA , 


( 6 ) 


where A is the surface area. Performing the substitution and integrating results in the energy 
expression: 
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«t; doW <.V> j • „4 sb /„2 m 2\ 2 , 

KB "n I T-^7 ?j “»■ 


The static solution exists at a minimum energy: 


( 7 ) 


3U . n 

WZ ' 0 ’ 

nm 


This results in all but equal to zero. The displacement function is therefore: 


aT' (l+v)a 2 b 2 

"<’•» ‘ ‘TiTTTT cos ? C0! * 


( 8 ) 


(9) 


Solving for moments by differentiating W and using equation (3) results in the equations for 
element moments: 


M„ = aT 


cos 2^ 

D 


M.. = aT 


^ [’ • ffi] cos * 

; cosi f • 


aT* D(l-\> )ab 


xy 


V -- 7— sin SL sin 22. 
a 2 + b 2 a b 


(10) 


(ID 


( 12 ) 


D. Results 

Figure 2 compares the element forces given by the above equation and the NASTRAN results. 
Figure 3 compares the normal displacements. The maximum errors for displacements, constraint 
forces, element forces and element stresses are listed in Table 1. 
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Figure 1. Simply-supported rectangular plate with a thermal gradient. 
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Table 1. Comparison of analytical and NASTRAN results. 


CATEGORY 

MAXIMUM 

ANALYTICAL 

MAXIMUM 

DIFFERENCE 

PER CENT 
ERROR 

Displacement 

6.2898 x 10' 1 

-1.5464 x 10' 3 

-0.25 

Constraint Force 

150.0 

-.9594 

-0.65 

Element Mom. , M 

X 

1.4770 x 10 2 

-1.1767 

-0.80 

Element Stress 

7.764618 x 10 3 

-90.33275 

-1.16 
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Figure 3. Displacement at 




RIGID FORMAT No. 1 (APP HEAT), Heat Conduction Analysis 
Linear Steady State Heat Conduction Through a Washer 
Using Solid Elements (1-12-1) 

Linear Steady State Heat Conduction Through a Washer 
Using Ring Elements (1-12-2) 

A. Description 

This problem illustrates the capability of NASTRAN to solve heat conduction problems. The 
washer, shown in Figure 1, has a radial heat conduction with the temperature specified at the out- 
side and a film heat transfer condition at the inner edge. Due to symmetry about the axis and the 
assumption of negligible axial gradients, the temperature depends only upon the radius. 

B. Input 

The first NASTRAN model is shown in Figure 2. The solid elements (HEXA1 , HEXA2, WEDGE and 
TETRA) and boundary condition element (HBDY, type AREA4) are used. The conductivity of the material 
is specified on a MAT4 card. Temperatures are specified at the outer boundary with SPC cards. 

Punched temperature output is placed on TEMP bulk data cards suitable for static analysis. 

Another variation of the problem is shown in Figure 3. Solid of revolution elements (TRIARG 
and TRAPRG) and boundary condition element (HBDY, type REV) are used. The conductivity of the 
material and the convective film coefficient are specified on a MAT4 card The CHBDY card references 
a scalar point at which the ambient temperature is specified using an SPC card. An SPC1 card is 
used to constrain the temperature to zero degrees at gridpoints on the outer surface. 


C. Theory 

The mathematical theory for the continuum is simple, and can be solved in closed form. The 
differential equation is 


1 

r 


JL (rk— ) = 0 
Br [rK dr } 


( 1 ) 


The boundary conditions are 

and 


'Hf = H(U a - u ) at r = r ] 
U = 0 at r = r 2 


( 2 ) 

(3) 
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The solution is 


U(r) 


HU a 

( k/ r-j > + H Mr z /r } l *- p(r 2 /r) 


= 288.516 *.n(2/r) 


D. Resul ts 

A comparison with the NASTRAN results is shown in Table 1. 


Table 1. Comparison of Theoretical and NASTRAN Temperatures for Heat Conduction in a Washer. 


r(radius) 

Theoretical 

Temperatures 

NASTRAN Temperatures 
(Sol ids)* 



NASTRAN Temperatures 
(Rinqs)* 

1.0 

199.984 

202.396 

199.932 

1.1 

172.486 

173.904 

172.448 

1.2 

147.381 

148.833 

147.355 

1.3 

124.288 

124.783 

124.269 

1.4 

102.906 

102.852 

102.894 

1.5 

83.001 

82.913 

82.992 

1.6 . 

64.380 

64.306 

64.375 

1.7 

46.889 

46.832 

46.886 

1.8 

30.398 

30.356 

30.397 

1.9 

14.799 

14.773 

14.798 

2.0 

0.000 

0.000 

0.000 


*These are the average temperatures at a radius. 
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Film heat transfer, 



Figure 1. Washer Analyzed in Heat Conduction Demonstration Problem 
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T TRIARG elements 
Q TRAPRG elements 
U = 488.5 at left end 

a 

U_ = 0.0 at right end 

d 


Figure 3. Section of a pipe, modeled with ring elements 
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RIGID FORMAT No, 1, Static Analysis 

Thermal and Pressure Loads on a Long Pipe Using Linear Isoparametric El^uents (1-13-1) 

Thermal and Pressure Loads on a Long Pipe Using Quadratic Isoparametric Elements (1-13-2) 

Thermal and Pressure Loads on a Long Pipe Using Cubic Isoparametric Elements (1-13-3) 

A. Description 

These problems demonstrate the use of the linear, quadratic and cubic isoparametric solid 
elements, IHEX1 , IHEX2 and IHEX3, respectively. A long pipe, assumed to be m a state of plane 
strain, is subjected to an internal pressure and a thermal gradient in the radial direction The 
structure modeled is shown in Figure 1. The finite element NASTRAN models for each of the elements 
are shown in Figures 2, 3 and 4. 

B. Input 

1. Parameters: 
r inner = a = 4 in. 
r outer = b = 5 in. 

E = 30. xlO 6 psi 
v = 0.3 

a = 1.428 x 10' 5 

p = 7.535 x 10‘ 4 1 - - sec - 
in 4 

p = 10 psi 
T = 100. 0°F 
T q = 0.0°F 

2. Boundary Conditions: 

u e = 0 at all points on the right side 

Ug = 0 at all points on the left side 

u z = 0 at all points on the bottom surface 

u z = 0 at all points on the top surface 


(radius to the inner surface) 
(radius to the outer surface) 
(Young's Modulus) 

(Poisson's Ratio) 

(thermal expansion coefficient) 
(mass density) 

(inner surface pressure) 

(inner surface temperature) 
(outer surface temperature) 
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3. Loads: 


Subcase 1 , 

p = 10 psi (internal pressure) 
Subcase 2, 


T 


r 


<yy 


ln(£) 


100 

ln(l .25) 


ln(£) 


where r is any radius. 


C. Theory 

1. Subcase 1 

The normal stresses due to the pressure load (Reference 24) are obtained by 

a 2 b 2 £ + pa 2 
" (b 2 -a 2 ) r 2 (b 2 -a 2 ) 

a 2 b 2 £ + pa 2 
(b 2 -a 2 ) r 2 (b 2 -a 2 ) 

(b 2 -a 2 ) 

where r is the radius and all shearing stresses are zero. 

The displacement in the radial direction is 

u = (1-2v)(l+v) pa 2 + (1+v) 1 pa 2 b 2 

E (b 2 -a 2 ) E r (b 2 -a 2 ) ’ 

and all other displacements are zero. 


= 


and 


2. Subcase 2 

The stresses in the radial and tangential directions due to the thermal load 
(Reference 24) are given by 


and 


ctET, 


2(l-v)ln(f) 


aET i 


2(l-v)ln(f) 


1 n(~) 5 

r (b 2 -a 2 ) 


(1 


1 - In 0 s— 

r (b 2 -a 2 ) 


■7> '"<?>] ■ 

<’*7> '"<;>] 
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The stress in the axial direction is obtained via the procedure contained in the 
reference as 


aET i 

20-v)ln(|) 


v - 


2a 2 v 
( b 2 -a 2 ) 


1"(?) - 2 1 


»<?>] 


All shearing stresses are zero. 

The displacement in the radial direction is 


u 


r 



t i 





a 2 b 2 

2(b 2 -a 2 ) 



+ i 


ln(£) 


+ 1 + (l-2v) 



D. Results 

Representative displacements and stresses for the finite element results com- 
pared to theoretical predictions are plotted in Figures 5 and 6. Note that five 
IHEX1 elements were used along the radial thickness whereas one element was used for 
each of the IHEX2 and IHEX3 cases. Two values for the stress occur at the boundary 
of two adjacent IHEX1 elements resulting in a sawtooth pattern. 
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Figure 2. Model of section using forty IHEX1 elements. 
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26 



Figure 4. Model of section using one IHEX3 element. 
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(a) Radial deflections, pressure load. 



(b) Radial deflections, thermal load. 


Figure 5. Deflection comparisons. 
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RIGID FORMAT No. 1, Static Analysis 
.Static Analysis of a Beam Using General .Elements (1-14-1) 


A. Description 

This problem demonstrates the use of general GENEL elements having various types of input 
format in the static analysis of a cantilever beam subjected to tension and bending. The beam 
consists of five GENEL elements and one BAR element as shown in Figure 1 
The GENEL elements are constructed as follows: 


GENEL Element 

Approach 

Matrix Size 

<V 

[S] 

1 

Flexibil ity 

3 

No 

No 

2 

Stiffness 

6 

No 

No 

3 

Stiffness 

3 

Yes 

Yes 

4 

Stiffness 

3 

Yes 

No 

5 

Flexibil ity 

3 

Yes 

No 


B. Input 

1 . Parameters 

£ = 6.0 m (length) 

E = 6.0 N/m 2 (modulus of elasticity) 

V = 0.3 (Poisson’s ratio) 

A = 1.0 m 2 (cross-sectional area) 

I = .083 m M (bending moment of inertia) 

F x = 1.0 N (axial load) 

P = 1.0 N (transverse load) 

C. Theory 

The stiffness matrix for the BAR element in its general form is given in section 8 of the 
NASTRAN Programmer’s Manual. 

Define [Z] as the matrix of deflection influence coefficients (flexibility matrix) whose terms 
are {u^} when {u^} is rigidly constrained, 

[K] as the stiffness matrix, 

[S] as a rigid body matrix whose terms are {^1 due to unit motions of {u d } , when all {^1 = 0, 
{f^ } as the vector of forces applied to the element at {u^ } , 
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and { f d > as the vector of forces applied to the element at {u d >. They are assumed to be statically 
related to the {f . } forces, i.e., they constitute a nonredundant set of reactions for the element. 
If transverse shear is neglected and the beam is confined to motion in the X-Y plane, then 

{f,l = [K] (u.) , 


where 


<V = 


{Ul } = 


6y 

9z 



LUl 

0 

o 


VO 

1 

0 

i 

o 

[K] = 

0 

12EI 

£3 

6EI 

£ 2 

= 

0 

6 

3 


0 

6EI 

£ 2 

4EI 
£ J 


1 

o 

3 

PO 

1 


and 





r I 
6 

0 

0 



[F] = 

[K]' 1 

0 

2 

3 

-1 

5 




_ 0 

-1 

2 



1 

0 

> 

C 

<< 


1 

0 

0 

[S] = 

0 

1 

Au x 

= 

0 

1 

-1 


_0 

0 

1 


0 

0 

1 _ 


where Au = u d - , i.e., the difference between the dependent displacement degree of freedom 

(u d > and the independent displacement degree of freedom 

D. Results 

The theoretical maximum deflection of the cantilever beam subjected to tension and bending 
(for the input values) are 

6x = ~~ = 1.0m (tension) 

and 
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These results are obtained by NASTRAN. 


n 


144. Q m (bending) 
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GENEL elements 1 thru 5 
R0D element 6 


Figure 1. NASTRAN General Element model. 
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RIOID FORMAT No. 1, Static Analysis 

Axisymmetric Cylinderical Thick Shell Subjected to Asymmetric Pressure Loadinq 

(1-15-1) 


A. Description 

This problem demonstrates the use- of elements TRAPAX and TRIAAX in the analysis of asymmetri- 
cally loaded solids of revolution The structure, illustrated in Fiqure 1, consists of a circular 
cylindrical shell loaded with a uniform external pressure over a small square area. 

The cylindrical shell wall is assumed to be simply suoported, i.e., the radial and circumferen- 
tial deflections and the bendino moments are zero at the ends. 

The upper half of the structure is modeled as shown in Figure 2. Trapezoidal elements havinq 
small and large dimensions, are used in the vicinity of the load and away from the load, respectively. 
A transition area, between the two trapezoidal configurations, is modeled with trianqular elements 
to illustrate their use. 

The loads and deflections, not required to be axi symmetric , are exoanded in Fourier series 
with respect to the azimuthal coordinate. Hue to the one plane of symmetry of this problem with 
resoect to the $ = 0 Diane, the deflections are represented by a cosine series selected bv the 
AXISY.M Case Control card. The highest harmonic used, in, is defined on the A*IC Pulk Data card. 

The oressure load is defined usinn PRESAX bulk Data cards. 

R. Input 

1. Parameters: 


r a 

= 

15 in. 

(Average radius) 

t 

= 

1 in. 

(Thickness) 

i 

= 

45 in. 

(Length) 

?c 

= 

3.75 in. 

(Load Lenoth) 

3 

= 

0.125 radians 

(Load Arc (3 = c/r )) 
a 

E 

= 

66666.7 psi 

(Modulus of Elasticity) 

V 

= 

0.3 

(Poisson's Ratio) 

n 

= 

10 

(Harmonics) 

Loads ' 



D 

= 

7.11111 psi 

(Pressure) 

A 

= 

14.063 in 2 

(Area of Load (A = 4c 
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3. Supports: 


Simply supported at the ends* u r = 0, u^ = 0 
Symmetry at the midplane: u z = 0 


C. Theory 


Theoretical results for this problem are taken from Reference 30, d. 568. The following 
theoretical values occur at the center of the load (z = $ = 0): 

u r = Er~ = ^.0272 in. (Radial Deflection (inward)) 

a 

= 0.1324 pA = 13. 24 in-lb/in (Circumferential Bendinq Moment) 

M = 0.1057 pA = 10.57 in-lb/in (Lonqi tudinal Bendinq Moment) 

= -2.6125 = -17.42 lb/in (Ci rcunferential Membrane Force) 

a 

F z = -2.320 = -15 . A 7 Ih/in (Lonmtudinal Membrane Force) 

a 

Theoretical stresses on the inside and outside walls at this point (z = p 4> = 0) are 

computed as follows: 


a 


z 



+ 


6M. 

4.2 


47.95 psi (Inside Wall Longitudinal Stress) 
-78.89 psi (Outside Wall Longitudinal Stress) 

62.02 psi (Inside Wall Circumferential Stress) 
-96.86 psi (Outside Wall Ci rcumferential Stress) 


D. Resul ts 

Figure 3 shows the NASTRAN radial deflection at the center of the load as a function of the 
number of harmonics selected for the solution. As can be seen, the solution is near convergence 
with ten harmonics. 

Figure 4 shows stresses, and o^, through the shell wall, at the center of the load. Ten 
harmonics shows very good convergence to nearly the theoretical values computed above. However, 
seven harmonics would result m relatively poor convergence even through Figure 3 indicates the 
displacement was close to convergence. Thus, displacement convergence alone may be an invalid 
indicator of an adequate solution. 
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Figure 1. Cylindrical shell loaded by a uniformly distributed load 
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Simply supported 



Symmetric 








Stress, psi Stress, psi 




(b) Circumferential stress. 


Figure 4. Stresses at center of load. 
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RIGID FORMAT No. 1, Static Analysis 
Fully Stressed Design of a Plate with a Reinforced Hole (1-16-1) 


A. Description 

A flat plate with a reinforced hole in the center is optimized for stresses due to a uniform 
end load. Restrictions on the minimum thickness are maintained. The plate is shown in Figure 1 
and the finite element idealization is illustrated in Figure 2. This problem has been investigated 
by G. G. Pope (Reference 21). 

Due to symmetry, only one quadrant is modeled. Due to the membrane load all rotations and 
normal displacements are constrained. The QDMEM and TRMEM elements are used for the plate and R0D 
elements for the reinforcement around the hole. 

The problem demonstrates several features unique to fully stressed design capability in NASTRAN. 
These features are: 

1. Elements with no limits on the range of the property change, i.e., the R0D has no 
PLIMIT data. 

2. Elements with a lower limit on the property optimization card. All membrane elements are 
required to have a resultant thickness which must not be less than a minimum thickness. 

This minimum is determined from the thickness obtained when the plate without a hole is 
subjected to an end load at a prescribed stress limit. 

3. Elements whose stress is not inspected but being in an area of nearly uniform stress have 

their properties changed due to another element's stress. Element 7 has no stress request 
but does have the same property identification number as element 17. This type of optimiza- 
tion can save computer time at the expense of a design that may not be truely optimized. 

4. A property whose value depends on the maximum stress of two elements. Elements 5 and 

15 have the same property card. This option may be necessary if insufficient core is 
allocated. 

5. Temperature dependent stress limits for material 3 

6. Using one stress limit only. The membrane elements use the maximum principle shear 

only. This is 1/2 the major principle stress allowed. This stress limit was chosen 

to better model the octahedral limit in Reference 21. 

The rod elements use only the tension and compression stress appropriate to the given 
property, namely area. 

7. An additional load case that was not included in the fully stressed design because a 
stress request was not made. The second subcase may be considered a displacement 
verification of this load case. 
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B. Input 


1. Parameters 
l = 30.0 in 
w = 20.0 in 
d = 10.0 in 
t * 3.348 in 


(total length) 

(total width) 

(hole diameter) 

(initial plate thickness) 

(initial rod cross sectional area) 


A q = 1.674 in 
E = 30.xl0 6 psi (modulus of elasticity) 
v = 0.3 


1.0 in 


(Poisson's ratio) 

(lower limit for plate thickness corresponding to a 25.0xl0 3 maximum 
principle stress) 


2. Boundary conditions: 

on y = 0 plane, u^ = 0 (symmetry) 
on x = 0 plane, u x = 0 (symmetry) 

all points u = 0 = 0 = 6 = 0 (permanent constraints) 

Z X jr Z 

3. Loads* 

First subcase: uniform load, F 1Q = 25.0xl0 3 lb/in 

Second subcase: at grid points 69 and 79, F^ = -1000.0 lb 

at grid point 78, F^ = -2000.0 lb 


(contact load on rim of hole - displacement check only) 


C. Theory 

The theoretical approach developed for the property optimization technique in NASTRAN is 
contained in the NASTRAN Theoretical Manual, Section 4.4. This technique is a fully stressed 
design approach. A mathematical programming technique is used in reference 21 from which the 
example problem was taken. 

The two techniques might be expected to give similar results when the same model is used. 
However, reference 21 employs elements which allow varying properties and stresses while NASTRAN 
elements allow only constant properties and constant stresses. Somewhat different geometry is 
used in the NASTRAN model, i.e., the use of quadrilateral elements for illustration. Additional 
features of the NASTRAN model are discussed in items 3, 4 and 5 of Part A. 


1.16-2 (3/1/76) 



D. Results 


The optimization process in this problem is terminated at 5 iterations. The initial weight to 
final weight ratio is 2.70 compared to Pope’s results of 2.63. Tables 1 and 2 show the optimized 
nondimensional properties of the elements around the arch. Note that the results from reference 
21 are averaged to provide an equivalent constant property element for comparison. 


Table 1. Optimized Nondimensional Thickness Comparisons. 


Element 

Original 

t/t e 

Reference 21 
Average 
t/t e 

NASTRAN 

37 

hem 

1.24 

1.00 

38 

1 

1.00 

1.04 

39 

3.348 

1.00 

1.00 

46 

3.348 

2.10 

1.14 

47 

3.348 

1.34 

2.00 

57 

3.348 

3.32 

1 . 34 

59 

3.348 

3.19 

4.40 

67 

3.348 

4.58 

5.47 

68 

3.348 

3.26 

1.00 

69 

3.348 

4.52 

5.49 


Table 2. Optimized Nondimensional Area Comparisons. 


Element 

Original 

A/dt e 

Reference 21 
Average 
A/dt e 

NASTRAN 

A/dt e 

■HP 

.1674 

.0249 

.00716 


.1674 

.0238 

0.0 effective 

mm 

.1674 

.0636 

.05019 

104 

.1674 

.1880 

.1839 

105 

.1674 

.3540 

.3287 
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Figure 1. Plate with reinforced hole. 
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Figure 2. Finite element model. 
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RIGID FORMAT No. 2, Inertia Relief Analysis 

Inertia Relief Analysis of a Circular Ring Under Concentrated and Centrifugal Loads (2-1-1) 

A. Description 

This problem illustrates the use of inertia relief analysis to solve a free-body problem. In 
inertia relief the structure is under constant acceleration due to the applied loads; the reac- 
tions to the applied load are due to the masses of the structure. Ficticious, nonredundant , sup- 
port points must be provided to define a reference system attached to the body. The displacements 
of the body are measured relative to the supported coordinates. 

The basic problem is illustrated in Figure 1. The structure consists of a spinning ring with 
a constant radial load applied to one point. The rotational velocity creates centrifugal loads 
and the point load causes inertia reactions. The actual dynamic motion of the whole structure 
is a cyclic motion of the center point coinciding with the rotation of the ring. The displace- 
ments measured by the inertia relief analysis, however, will be the static motion relative to the 
support point displacements. 

The displacements are defined in a cylindrical coordinate system (u-j = u r , u^ = u^, u^ = u z ). 
The elements used are BAR elements with a large cross-sectional area to minimize axial deforma- 
tions. The BARs were offset a uniform radial distance from the grid points to demonstrate the 
offset option of the BAR element. 

B, Input 

1 . Parameters : 

R = 10.0 (radius at end of BAR elements) 

R^ - 11.0 (Radius at grid points) 

I = 10.0 (Bending inertia) 

p = 0.5 (Mass density) 

E = 1000. (Modulus of elasticity) 

A = 1000. (Cross-sectional area) 

2. Loads: 

P r,13 " 100 

f = 1.59 cps (Rotational velocity, u = 1.0 radians per second) 
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3. Supports* 

a) The u direction is supported to restrict vertical translation. 

r » i 

b) The u^ i and u^ ^ directions are supported to restrict rotation and horizontal 
translation. 


4. Grid Point Weight Generator Input: 

Weight and moment of inertia are defined relative to point 19. 


C. Answers 


1. The Element Forces and Moments may be solved by the following analysis, as explained in 
Reference 7, Chapter 12. 

a) Using symmetry the structure may be defined by the free-body diagram in figure 2. 
The static equilibrium equations at any angle are: 

A = A q coscj) + y<j> sin<f> (Axial Force) 

V = A q sin<j> + y<j> cos (|> (Shear) 

M = M q + r[y(l - cos<f> - <f> si ncj> ) + A q (1 - cos<|>)] (Bending Moment) 

b) Using energy and Castigli ano‘s Theorem: 


U 

6U 

o 


R % r m 2 , , 

2eT / M d<)> 

0 

0 


6U 

6A o 


0 


These are the deflections at the bottom which are fixed. The resulting two equa- 
tions are used in step c. 

c) Solving the equations in (b) gives the redundant forces- 


M 


o 



FR 

4 
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d) Adding a dummy load and solving the problem with the above boundary conditions gives 
the displacement due to the point load: 

6 - fr3 r* 2 n 

6 f ~ an ( §- - 1} 

e) The axial stress and radial displacement due to the centrifugal load is: 

a u = pR 2 to 2 = 5.0 x 10 2 


9 2 

2pR w 


= 1.0 


f) The total result of summing the two loads is: 



THEORY 

NASTRAN 

Displacement u^ ^ 

1.75 

1.734 

Moment BAR #1 , end A 

-79.5 

-80.48 

Moment BAR #12, end B 

-238.5 

-236.0 


g) The structural mass characteristics as calculated by the grid point weight 
generator are: 


Theoreti cal 


X c q =11 .0 from point 19 
Mass = if x 10 4 = 3.14159 x 10 4 


I = I = t X 10 6 = 1.5708 x 10 6 

xx yy 2 

l zz = it x 10 6 = 3.14159 x 10 6 


NASTRAN 

11.0 

3.1326 x 10 4 
1.5663 x 10 6 
3.1326 x 10 6 


(Inertias are about center of gravity) 

NASTRAN gives slightly different answers due to the polygonal shape of the finite element model. 
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Figure 2. 


Free body diagram of loads in bending ring. 
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RIGID FORMAT NO. 2, Inertia Relief Analysis 
Windmill Panel Sections for Automated Multi-stage Substructuring , Run 1, (2-2-1) 

Windmill Panel Sections for Automated Multi-stage Substructuring, Run 2, (2-2-2) 

Windmill Panel Sections for Automated Multi-stage Substructuring, Run 3, (2-2-3) 

Windmill Panel Sections for Automated Multi-stage Substructuring, Run 4, (2-2-4) 

Windmill Panel Sections for Automated Multi-stage Substructuring, Run 5, (2-2-5) 

Windmill Panel Sections for Automated Multi-stage Substructuring, Run 6, (2-2-6) 

Windmill Panel Sections for Automated Multi-stage Substructuring, Run 7, (2-2-7) 

A. Descri ption 

This problem illustrates the fully automated multi-stage substructuring capability of NASTRAN. 
The single structure model for the Windmill panel problem is shown in Figure 1. Indicated in this 
figure are the three basic substructures used for the analysis. As can be seen, the entire structure 
can be composed of only these three components, thus taking advantage of symmetry. The detailed 
idealizations for the three basic substructures are shown in Figures 2 and 3. These figures show 

the three separate basic coordinate systems and the local coordinate systems for each of the three 

basic substructures created. 

Of the total of seven runs involved, three Phase 1 runs are made, one for each basic substruc- 
ture, using Rigid Format 2 in order to generate mass matrices. The combination and reduction to 
the final model is accomplished in seven distinct Phase 2 steps, plus eight equivalence operations. 
The sequence of combination steps taken is illustrated in Fiaures 4a and 4b. Figure 5 details the 

points retained in the "analysis set" following the Phase 2 Guyan reduction. A static solution, 

Rigid Format 1, is obtained for each of the three load cases specified. Run 4 produces actual 
plot output. Runs 5 and 6 demonstrate the Phase 3 data recovery for two of the basic substructures . 

A seventh run is made to extract normal modes using Rigid Format 3 for the same reduced struc- 
ture shown in Figure 5. 

B, Input 

1. Parameters. 

r Q = 50.0 in (outer radius) 

r^ = 10.0 in (inner radius) 

t = 0.1 in (plate thickness) 

E = 10 x 1,0^ psi (modulus of elasticity) 

v = 0.25 (Poisson's ratio) 
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2. Boundary Conditions: 

All points u - e = e = 0 = 0 (permanent constraint) 

^ a Jr Z 

u x = 0 at HUB grid points 13, 19, 37, 43 
Uy = 0 at HUB grid points 1, 7, 25, 31 

3. Loads: 

First Subcase: centrifugal force due to unit angular velocity 

Second Subcase: unsymmetric load - right panel in tension, bottom panel in 

compression, F = 100 uniformly distributed over each loaded edge 
Third Subcase: F = 1.0 applied at HUB grid point 4 inward radially 

4. Substructuring Parameters: 

S0F ( 1 ) = S0FO,95O $ CDC 
S0F ( 1 ) = FT18,950 $ IBM 
S0F(1 ) = INPT,950 $ UNIVAC 
PASSWORD = DEM0 
0PTI0NS = K, M, P 


C. Theory 

This problem is designed to illustrate the use of automated multi-stage substructuring. No 
closed form solution is available. Results are compared with non-substructured NASTRAN solutions. 

D. Results 

The solutions of the final reduced structure using both Rigid Format 1 and Rigid Format 3 are 
in excellent agreement with the non-substructured solutions. Displacements at selected points and 
eigenvalues are compared in Table 1. The values presented were obtained from executions on IBM 
equipment. Values obtained from CDC and UNIVAC are of the same order of magnitude with slight 
differences attributable to round-off of very small numbers. 
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Figure 3. Windmill section substructures. 
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RIGID FORMAT No. 3, Real Eigenvalue Analysis 

Vibration of a 10x20 Plate (3-1-1) 
Vibration of a 20x40 Plate (3-1-2) 
Vibration of a 10x20 Plate (INPUT, 3-1-3) 
Vibration of a 20x40 Plate (INPUT, 3-1-4) 


A. Description 

This problem demonstrates the solution for natural frequencies of a large-order problem. The 
structural problem consists of a square plate with hinged supports on all boundaries. The 10x20 
model, as shown in Figure 1, uses one half of the structure and symmetric boundary constraints on 
the mid- line in order to reduce the order of the problem and the bandwidth by one-half. The 20x40 
model is the same dimensions with a finer mesh. Both configurations are duplicated via the INPUT 
module to generate the QUAD1 elements. 

Because only the bending modes are desired, the in-plane deflections and rotations normal to 
the plane are constrained. The inverse power method of eigenvalue extraction is selected and both 
structural mass density and non-structural mass-per-area are used to define the mass matrix. 

Table 1 lists the NASTRAN and theoretical natural frequencies as defined in Reference 8. 
Figures 2 and 3 are comparisons of the first two mode shapes. The modal masses for these modes 
»are equal to one-fourth the total mass or m.. = 10302.2 . 

An undeformed structure plot is executed without plot elements. This is expensive on most 
plotters since all four sides of each quadrilateral are drawn. Plot elements are used to draw 
an edge only once and to draw selected coordinate lines (every second or fourth line depending on 
the model used) for the deformed plots of each eigenvector. 

B. Input 

1 . Parameters : 


l 

= w = 20.0 

(Length and width) 

I 

_ 1 
' 12 

(Moment of inertia) 

t 

= 1.0 

(Thickness) 

E 

= 3 x 10 7 

(Modulus of elasticity) 

V 

= 0.30 

(Poisson's ratio) 

P 

= 206.0439 

(Mass density) 
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2. Boundary constraints: 


along x = 0, 0 y = 0 

along y = 0, u z = 0 y * 0 

along x = 10, u z = 0 x = 0 

along y = 20, u z = 0 y = 0 

3. Eigenvalue extraction data: 

Method: Inverse power 

Region of interest: .89 < f < 1.0 

Number of desired roots: 3 

Number of estimated roots: 1 


Symmetric Boundary 


Hinged Supports 


Table 1. Natural Frequency Comparisons, cps. 


Mode 

No. 

Theoretical 

N AST RAN 
10x20 

N AST RAN 
20x40 

i 

.9069 

.9056 

.9066 

2 

2.2672 

2.2634 


3 

4.5345 

4.5329 
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RIGID FORMAT No. 3, Real Eigenvalue Analysis 
Vibration of a Compressible Gas in a Rigid Spherical Tank (3-2-1) 


A. Description 

This problem demonstrates a compressible gas in a rigid spherical container. In NASTRAN a 
rigid boundary is the default for the fluid and, as such, no elements or boundary lists are 
necessary to model the container. 

Aside from the NASTRAN bulk data cards currently implemented, this problem demonstrates the 
use of the hydroelastic data cards: AXIF, CFLUID2 , CFLUID3, and RINGFL. 

i 

The lowest mode frequencies and their mode shapes for n = 0, 1 and 2 are analyzed where n is 
the Fourier harmonic number. Only the cosine series is analyzed. 

B. Model 

1 . Parameters 

R = 10.0 m (Radius of sphere) 

p = 1.0 x 10”^ Kg/rn^ (Mass density of fluid) 

B = 1.0 x 10^ Newton/m^ (Bulk modulus of fluid) 

2. Figure 1 and 2 show the finite element model. The last 3 digits of the RINGFL identifica- 
tion number correspond approximately to the angle, 0, from the polar axis along a 
meridian. 

C. Theory 

From Reference 18, the pressure in the cylinder is proportional to a series of functions: 


Q 


n,m 


J m+i < X > n 

— i P n (cos 6) cos ncp , 

/T m 


n < m 

m = 0,1 ,2 


( 1 ) 


where: 


(3 Pressure coefficient for each mode 
n ,m 

w mk 

X Nondimensional radius = -— r 

a 

w . Natural frequency for the kth mode number and mth radial number in radians 
per second 

Bessel function of the first kind 
m+i 
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r 


radius 


a 



speed of sound in the gas 


associated Legendre functions 

0 meridinal angle 

<f) circumferential angle 

n harmonic number 

m number of radial node lines 


The solution for X and hence is found by the use of the boundary condition that the flow 
through the container is zero. 


W x) 

L A 


= 0.0 


( 2 ) 


/r=R 


where R is the outer radius . i ; 

! 

This results in zero frequency for the first root. Multiple roots for other modes can be 
seen in Table 1. The finite element model assumes different pressure distributions in the two 
angular directions which causes the difference in frequencies. 


D. Resul ts 

Table 1 and Figure 3 summarize the NASTRAN and analytic results for the lowest nonzero root 
in each harmonic. Table 1 lists the theoretical natural frequencies, the NASTRAN frequencies, the 
percent error in frequency, and the maximum percent error in pressure at the wall as compared to 
the maximum value. Figure 3 shows the distribution of the harmonic pressure at the wall versus 
the meridinal angle. The theoretical pressure distributions correspond to the Legendre functions 
P Q (cos 0), P Q (cos 6), and P Q (cos 6) which are proportional to cos 0, sin e, and sin 0 respec- 
tively. 
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Table 1. Comparison of NASTRAN and analytical results. 


— 

Harmonic 

Natural Frequency (Hertz) 

Pressure 

Analytical 

NASTRAN 

% Error 


0 

33.1279 

33.2383 

0.33 

1.19 

1 

33.1279 

33.2060 

0.24 

0.47 

2 

53.1915 

53.3352 

0.27 

0.91 
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Figure 2. Gas filled rigid spherical tank model. 
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RIGID FORMAT No. 3, Real Eigenvalue Analysis 
Vibration of a Liquid in a Half-Filled Rigid Sphere (3-3-1) 


A. Description 

The model is similar to Demonstration Problem No. 3-2 except that a hemispherical fluid model 
with a free surface is analyzed. Additional cards demonstrated are the free surface list (FSLIST) 
and free surface points (FREEPT). The effective gravity for the fluid is found or the AXIF card. 
The fluid is considered incompressible. 

The lowest three eigenvalues and eigenvectors for the cosine and sine series of n = 1 are 
analyzed, where n is the harmonic order. 

B. Input 

1 . Parameters 

g = 10.0 ft/sec 2 (Gravity) 

R = 10.0 ft (Radius of hemisphere) 

p = 1.255014 lb-sec 2 /ft 4 (Fluid mass density) 

B = 00 (Bulk modulus of fluid, incompressible) 

2. Figure 1 shows the finite element model. 

C. Results 

Reference 17 gives the derivations and analytical results. In particular, the parameters 
used in the reference are: 

e = 0 (half-filled sphere) , ) 

> 0 ) 

2 \ 

X = (dimensionless eigenvalue) . ) 

Table 2 of Reference 17 lists the eigenvalues, X-| , and X^ for the first three modes. 
Figure 13 of Reference 17 shows the mode shapes. 

The analytic and NASTRAN results are compared in Table 1. The frequencies are listed and 
the resulting percentage errors are given. The maximum percent error of the surface displacement, 
relative to the largest displacement, is tabulated for each mode. 
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The free surface displacements may be obtained by the equation: 


where p is the pressure at the free surface recorded in the NASTRAN output. Note that, since an 
Eulerian reference frame is used, the pressure at the original (undisturbed) surface is equal to 
the gravity head produced by motions of the surface. Special FREEPT data cards could also have 
been used for output. Since the results are scaled for normalization anyway, the harmonic pres- 
sures may be used directly as displacements. 

Figure 2 is a graph of the shape of the free surface for each distinct root. Both analytic 
and NASTRAN results are scaled to unit maximum displacements. Because the cosine series and the 
sine senes produce identical eigenvalues, the resulting eigenvectors may be linear combinations 
of both series. In other words the points of maximum displacement will not necessarily occur at 
<f> = 0° or <f) = 90°. Since the results are scaled, however, and the results at <j> = 0 are propor- 
tional to the results at any other angle, the results at <f> = 0 were used. 


Table 1. Comparison of natural frequencies and free surface mode shapes from the reference and 

NASTRAN. 


Mode 

Natural 

Frequency (Hertz) 

Mode Shape 

Number 

Reference 

NASTRAN 

NASTRAN 
% Error 

Maximum 
% Error, e 

1 

0.1991 

0.1988 

-0.1 

e < 1 % 

2 

0.3678 | 

0.3691 

0.3 

e < 2.6 % 

3 

0.4684 

0.4766 

1.8 

e < 4 % 
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Figure 1. Rigid sphere half filled with a liquid. 
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RIGID FORMAT No. 3, Real Eigenvalue Analysis 
Acoustic Cavity Analysis (3-4-1) 


A. Description 

This problem illustrates the use of NASTRAN to determine the acoustic modes in a cavity 
containing both axisymmetric regions and evenly spaced radial slots. The solution is based 
on an analogy between pressure and di splacement , and between fluid particle acceleration and 
internal structural force described in the Theoretical Manual. 

B. Input 

The finite element model for the motor cavity of the Minuteman III, Stage III, is shown in 
Figure 1 As may be seen, it consists of six slots and a long, slender central cavity of 
irregular shape. The model consists of AXIF2 , AXIF3 , and AXIF4 finite elements in the central 
cavity, and SL0T3 and SL0T4 finite elements in the slotted region 

C. Results 

The vibration mode frequencies for harmonic n = 0 as determined with NASTRAN are shown in 
Table 1. Also shown are the vibration mode frequencies as determined with an acoustic model 
and reported in Reference 19. 
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Table 1. Natural frequencies for the third stage, Minuteman III, motor cavity. 


Mode 

Frequency, Hz 

NASTRAN 

Experi- 

mental 

1 

0.0 

0.0 

2 

90.1 

93.0 

3 

199.5 

200.0 

4 

310.4 

312.0 

5 

388.0 

388.0 

6 

449.1 

466.0 

7 

512.8 

518.0 
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RIGID FORMAT No. 3 (APP HEAT), Nonlinear Heat Conduction 
Nonlinear Heat Transfer in an Infinite Slab (3-5-1) 


A. Description 

This problem demonstrates NASTRAN's capability to solve nonlinear steady state heat conduction 
problems. The infinite slab is subjected to uniform heat addition per unit volume. There is no 
heat flux on one face and the other face is kept at zero degrees. The conductivity is temperature 
dependent. This is a one dimensional problem, since there is no temperature gradient parallel to 
the surfaces of the slab 


B. Input 

The NASTRAN model is shown in Figure 1 Linear elements BAR, C0NR0D, R0D and TUBE with areas 
of t r square units and boundary condition element HBDY (P0INT) are used. The heat addition is spec- 
ified on a QV0L card and is referenced in Case Control by a L0AD card. The area factor for the 
HBDY is given on the PHBDY card and heat flux is zero. The initial temperatures are qiven on a 
TEMPD card and referenced in Case Control by a TEMP (MATERIAL) card. The conductivity is specified 
on a MAT4 card and is made temperature dependent by the MATT4 card referencing table TABLEM3. The 
convergence parameter, the maximum number of iterations and an option to have the residual vector 
output are specified on PARAM cards. The temperature at the outer surface is specified by an SPC 
card. Temperature output is punched on TEMP bulk data cards for future, use in static analysis 


C. Theory 

The conductivity, k, is defined by 

k(T) = 1 + T/100 , 

where T is the temperature. 

The heat flow per area, q, is 

q(x) = -k ^ = -(1 + T/100) g . 

The heat input per volume, q y , affects the heat flow by the equation 


A convenient substitution of variables in Equations (1) and (2) is 

u = -/q(x)dx = (T + T 2 /200) . 

Differentiation and substitution for q in Equation (2) results in the second-order equation 
in u: 


( 1 ) 

( 2 ) 

(3) 
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( 4 ) 


A- a 

2 " q v * 
dx 


From the following boundary conditions 


and 


u = 0 at x = l 


&■» at ■ 


(5) 

( 6 ) 


the solution to Equation (4) is 


u = y U 2 - x 2 ) 


Therefore the solution for the temperature is 


T- = 100 [-1 + (l+q v U 2 - x 2 )/100) 2 ] 


(7) 


Since heat is flowing into the system, the positive temperature solution will occur. 


D. Results 

A comparison with NASTRAN results is shown in Table 1. 

Table 1. Comparison of theoretical and NASTRAN temperatures for nonlinear heat conduction in 
an infinite slab. 


Grid 

Point 

Theoretical 

Temperature 

NASTRAN 

Solution 

1 

73.20 

73.13 

2 

69.56 

69.53 

3 

58.11 

58.11 

4 

36.93 

36.93 

5 

0.00 

0.00 
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Grid 1 Flux = 0.0 
Grid 5 Tempe.ature = 0.0 


Figure 1. Slab modeled with linear elements 
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RIGID FORMAT No. 3, Approach Heat, 

Nonlinear Radiation and Conduction of a Cylinder (3-6-1) 


A. Description 

This problem illustrates the solution of a combined conduction and radiation heat transfer 
analysis. The model is a two-dimensional representation of a long cylinder subject to radiant 
heat from a distant source. The shell has internal radiation exchange, external radiation loss, 
and conduction around the circumference. 


B. Input 

The NASTRAN Model, shown in Figure 1, uses R0D elements to represent the circumferential heat 
flow and HBDY elements to represent the inside and outside surfaces. The radiation exchange 
factors for the inside of the cylinder are defined on the RADMTX data cards. The incominq vector 
flux is defined on the QVECT data card. The model parameters are: 


R = 2.0 ft 

t = .001 ft 

A = 20.306 ft 

e = a = 0.1 

q y = 425 BTU/(ft 2 -hr) 

k = 94.5 BTU/(hr-ft-°F) 

i-8 


(Radius of shell ) 

(Thickness) 

(Axial length) 

(Emissivity and absorptivity) 
(Source flux density) 
(Conductivity of shell) 


174 x 10”° BTU/(ft -hr-°R ) (Stefan-Bol tzmann radiation constant) 


C. Theory 

A closed-form solution to this problem is not available. However, the solution may be 
validated by checking the global net heat flow, the local net heat exchange, and the estimated 
average temperature. 

An estimate of the average temperature may be obtained from the equations: 


% 


aq v £R 


tt/2 

cos e de 

-tt/2 


2aJlRq v , 


( 1 ) 


and 


Q out = eaf4 (2 * Ril) 


(2) 
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where Q in is the total input from the source, Q Qut is the net flux radiated outward and T is 
the average absolute temperature. 

Since the net heat flow must be zero in a steady-state analysis, Equations (1) and (2) are 
equated to obtain: 


D. Results 

The resulting temperature distribution around the circumference of the shell is shown in 
Figure 2. The average value of temperature from the NASTRAN results shows 57.87° F. The esti- 
mated average temperature from Equation (3) above is 68°. The difference is due to the non-uniform 
radiation effects. 

A second check is provided by computing the global net heat flow error in the system. Summinq 
the net flow into each element gives a net heat flow error several orders of magnitude less than 
the total heat from the source. As a further check, the local net heat flow error at grid point 2 
was calculated by summing the contributions from the connected elements. The heat flow terms 
shown in Figure 3, as calculated by NASTRAN, were: 


q 2 

= 

59.420 

(Flow through R0D #2 (flux * area)) 

q 3 

= 

97.862 

(Flow through R0D #3 (flux • area)) 

Q r42 

= 

-133.564 

(Inside radiation flow into HBDY #42) 

Q r43 

= 

-85.352 

(Inside radiation flow into HBDY #43) 

Q r22 

= 

-305.418 

(Outside radiation into HBDY #22) 

Q r23 

= 

-257.930 

(Outside radiation into HBDY #23) 

Q v22 

= 

481.157 

(Vector flux input to HBDY #22) 

Q v23 

= 

381.848 

(Vector flux input to HBDY #23) 

The net flow 

error into grid point 2 is* 


Q 2 = 2 ( Q r22 + Q r23 + Q r42 + Q r43 + Q v22 + Q v23^ + Q 2 ' Q 3 

This error isless than 1% of the total heat flow input at the point. 


1.9 BTU (4) 
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Figure 1. Cross section of thin wall shell 
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ion heat transfer. 






Figure 3. 



O grid points 

rod elements 

— heat elements 


Illustration for heat exchange computation at a grid point. 
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RIGID FORMAT No. 4, Differential Stiffness Analysis 
Differential Stiffness Analysis for a Hanging Cable (4-1-1) 


A. Description 

Advanced versions of NASTRAN provide an iteration procedure for nonlinear differential stiff- 
ness (or geometric stiffness) solutions. As described in Section 7 of the NASTRAN Theoretical 
Manual, the internal loads are recalculated for each iteration. The chanqes in direction of these 
internal loads are used to correct the previous solution. External loads retain their oriqinal 
orientation, however, they do travel with the qrid ooint. 

A classical nonlinear qeonetric problem is that of a haminn cable which assumes the shape of 
a catenary when a uniform qravity load is aDDlied. As shown in Figure 1, the model is qiven a 
circular shape initially. The resulting displacements of the qrid points, when added to their 
onainal locations provide a close approximation to the catenary. 

B. Input 

The NASTRAN model consists of nine BAR elements connected to ten grid points evenly spaced 
on a quarter circle. The bending stiffness of the elements is a nominally small value necessary 
to provide a non-singular linear solution. 

The axial stiffness of the elements is a sufficiently large value to limit extensional dis- 
placements. The basic parameters are 

R = 10.0 ft. (initial radius), 

w = 1.288 lb/ ft (Weiqht per lenqth), 

and L = 5 tt. 

C. Theory 

With reference to the coordinate system illustrated on the next page, the basic differential 
equation, obtained from Reference 25 is 

& ■ 5(1 *<y ') 2 ) 1/2 , 0 ) 

where 

w is the weight per unit lenqth, 

H is the tension at x = 0, 

and y' = dy/dx is the slope of the resulting curve. 
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Dividing both sides of Equation (1) by the radical term and inteqratinq, results in the equation, 

slnh'V = 5 Jf +c 1 . (2) 

Since y' = 0 at x = 0, = 0 and 


y' = sinh (if) • 


(3) 


Integrating again and applyinq the known boundary condition y = 0 at x = 0, the equation for the 
shape is 


y- 


H 

W 


cosh (Sf)-l • 


w 


Since the length of the cable is known but the horizontal force, H, is unknown, the two may 
be related by integrating for the length, L, which is 


L 


- sinh 
w 



(5) 


where x Q is one-half the distance between supports. If w, x , and L are given, Equation (5) is 
solved for H (for x q = 10.0, w/H = .1719266) and Equation (4) is evaluated to obtain the actual 
shape. However, for purposes of comparison to the NASTRAN solution the location of several points 


fixed on the string are determined. For a given position, s, along the cable, the coordinates x 
and y would be 
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D. Results 


The following table compares theoretical results to those of NASTRAN. Deflections are 
measured from the initial shape at selected locations. 


Table 1. Comparison of NASTRAN Results to Theoretical Predictions. 


Grid 

Point 

s 

- Horizontal 

Uy - Vertical 

Theory 

NASTRAN 

Theory 

NASTRAN 

11 

13.962 

-.4856 

-.4739 

-.1119 

-.0408 

13 

10.472 

-.8043 

-.7666 

-.2286 

-.1269 

15 

6.981 

-.5175 

-.4612 

.0030 

.1470 

17 

3.491 

-.1110 

-.0877 

.5698 

.7973 

19 

.0 

.0 

.0 

.9338 | 

1.2167 
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(a) Initial quarter circle modeled. 



(b) Final catenary shape obtained. 


Figure 1. Hanging cable. 
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RIGID FORMAT No. 5, Buckling Analysis 
Synmetric .Buckling of a Cylinder (5-1-1) 


A. Descri ption 

This problem demonstrates the use of buckling analysis to extract the critical loads and the 
resulting displacements of a cylinder under axial loads. The Buckling Analysis rigid format solves 
the statics problem to obtain the internal loads in the elements. The internal loads define the 
differential stiffness matrix [K d ] which is proportional to the applied load. The load factors, 

X.. , which causes buckling are defined by the equation: 

[^•j [K d ] + [K]]{ Ui } = 0 

where [K] is the linear stiffness matrix. This equation is solved by the Real Eigenvalue Analysis 
methods for positive values of X. . The vectors {u^ are treated in the same manner as in real 
eigenvalue analysis. 

The problem is illustrated in Figure 1, it consists of a short, large radius cylinder under a 
purely axial compression load. A section of arc of 6 degrees is used to model the axisymmetric 
motions of the whole cylinder as shown in Figure 2. 

All three types of structure plots are requested: undeformed, static and modal deformed. The 
undeformed perspective plot is fully labeled for checkout of the problem. The modal orthographic 
plots specify a range of vectors {u..} which includes all roots. A longitudinal edge view of the 
model is also plotted for easy identification of mode shapes. 

B. Input 

1. Parameters' 


R = 

80 

(Radius) 

h = 

50 

(Height) 

E = 

1.0 x 10 4 

(Modulus of elasticity) 

v = 

0.0 

(Poissons rat o) 

t = 

2.5 

(Thickness) 

r b = 

1.30208 

(Bending inertia) 
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2 . 


Loads: 


p = 1.89745 x 10 3 /3° ARC 

3. Constraints: 

a) The center point (17) is constrained in u z . 

c) All points are constrained in u Q , 0 r# and 0 Z . 

d) The top and bottom edges are constrained in u r< 

4. Eigenvalue Extraction Data: 

a) Method: Unsymmetrical Determinant 

b) Region of Interest: .10 < A < 2.5 

c) Number of estimated roots =4 

d) Number of desired roots = 4 

e) Normalization: Maximum deflection 

C. Answer 

The solution to this problem is derived in Reference 9, p. 439. For axisymmetnc buckling, 
the number of half-waves which occur when the shell buckles at minimum load are* 


h /l 2 ( 1 - v 2 
71 V R 2 t 


where m is the closest integer to the right-hand values. 

The corresponding critical stress is: , 

Et 2 mV + Eh 2 
° cr 12h 2 (l-v 2 ) R 2 m 2 ir 2 

Using the values given, the lowest bulkling mode consists of a full sine wave. The NASTRAN results 
and the theoretical solutions for the critical load for each buckling mode are listed below- 


Number of 
Half Waves 

NASTRAN 

ANALYTICAL 

m 



1 

2.2889 

2.2978 

2 

.99424 

1.0 

3 

1.2744 

1.26402 

4 

2.U07C 

1 86420 
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Figure 1. Cylinder under axial load. 
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RIGID FORMAT No. 6, Piecewise Linear Analysis 
Piecewise Linear Analysis of a Cracked Plate (6-1-1) 


A. Description 


This problem illustrates elastic-plastic deformation of a thin plate uniaxially loaded across 
a crack at the center of the plate shown in Figure 1. The same problem was solved by J. L. Swedlow 
(Reference 10). 

Piecewise Linear Analysis involves loading the plate in increments and recalculating the 
material properties for each element as a function of the elemenl stresses for the last load 
increment. 

B. Input 

1. Parameters: 

L = 9.0 inch (Total length of plate) 

W = 6.0 inch (Total width of plate) 

w = 2.0 inch (Total width of crack) 

t = 1.0 inch (Thickness) 

E Q = 10.8 x 10 lb/in (Modulus of elasticity at zero strain) 

v Q = .3333 (Poisson's Ratio at zero strain) 

2. Loads: o is the applied load (Figure 1) 


Load Factor 

o lb/in 2 

Load Factor 

a lb/in 2 

1 

2,300 

14 

7,000 

2 

2,500 

15 

7,500 

3 

2,800 

16 

8,000 

4 

3,100 

17 

8,500 

5 

3,400 

18 

9,000 

6 

3,700 

19 

9,500 

7 

4,000 

20 

10,000 

8 

4,400 

21 

10,500 

9 

4,800 

22 

11,000 

10 

5,200 

23 

11,500 

11 

5,600 

24 

12,000 

12 

6,000 

25 

12,500 

13 

6,500 

26 

13,000 
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3. Constraints; 


a) All grid points are constrained in u , 6 , 6 , and 0 . 

z x y z 

b) Grid points along the Y-axis are constrained in the u x direction. 

c) Grid points along the X-axis from the crack tip (x = 1.0) to the edge (x = 3.0) are 
constrained in the u^ direction. 

C. Modeling Techniques 

The finite element model, shown in Figures 2 and 3, utilizes two planes of symmetry so only 
one quarter of the structure (the first quadrant) is modeled. All membrane elements use stress- 
dependent materials, duplicating the model in Reference 10. 

The octahedral stress used in the determination of the material properties was defined in 
Reference 10 as: 




°x ■ 


a a i + 
x y 


a 2 + 3a z 
y xy 


The octahedral strain was defined by. 


where 




£ P 


T limit 
a limi t 


9.716 x 10- 3 (T 0 /T„ m1t 

' (/3 ' 3 >°l»it 

= 11,500 lb/in 2 


1 ) 


1/0.3964 


T o ~ T limit 
T o > T limit 


NASTRAN uses an equivalent uniaxial stress-strain curve defined by 


a = 3//2 t q 

e = a/E + Jl Cp 
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This curve is shown in Figure 4. A complete discussion of the equations may be found in 
Reference 10. 

D. Answers 

Comparisons of analytical and calculated stresses in the elements along the axis of the crack 
are given in Figures 5 and 6. The analytical results, based on linear analysis, are compared with 
the calculated results at the end of the first load increment. 

Figures 7 through 9 show the displacement at the center of the crack and stresses for elements 
near the tip of the crack for all load factors. In the NASTRAN analysis, the octahedral stress is 
calculated for each load factor as a function of the current values of the stresses. In 
Reference 10 the current value of the octahedral stress is obtained by accumulating incremental 
values of the octahedral stress. This procedure results in a generally more flexible model as can 
be seen from the displacements in Figure 7. The resulting differences in calulated stresses are 
particularly noticeable at the higher load levels. 
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Figure 1. Cracked plate with uniaxial load. 
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Plane of Symmetry 



Figure 3. Triangular membrane element identification numbers near crack. 


- lb/in. 



6 . 1-7 ( 12 - 1 - 69 ) 


LO 



6.1-8 (6/1/72) 


Figure 5. Comparison of stress a along the X-axis, load factor 
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Figure 6. Comparison of stress a along the X-axis, load factor 



a 1 b/ i n . x 10 


Figure 7. Crack opening vs. load at center of crack (x 
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(a- 


5.5 



Figure 9. 


0y/a vs. a, 


four elements at tip of crack. 
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RIGID FORMAT No. 7, Complex Eigenvalue Analysis - Direct Formulation 

Complex Eigenvalue Analysis of a 500-Cell String (7-1-1) 

Complex Eigenvalue Analysis of a 500-Cell String ( INPUT, 7-1 -2) 

A. Description 

This problem demonstrates both the use of direct complex eigenvalue analysis and the various 
methods of supplying damping to the structure. The simulated model is a string under tension having 
uniform viscous and structural damping as shown in Figure 1. The stiffness due to tension is modeled 
with scalar springs, the mass is represented by scalar masses, and the viscous damping is provided by 
scalar dampers connected on one end to the points and fixed on the other end. The structural damp- 
ing is provided by the scalar springs and an overall damping factor, g^. The INPUT module is used 
to generate the scalar springs, dampers and masses. 

Complex Eigenvalue Analysis is used to solve the following general matrix equation: 

([M]p 2 + [B]p + [K] ) { u } = 0 


where 

p is the complex root 

[M] is thfe complex mass matrix 

[B] is the complex damping matrix for viscous damping 

[K] is the complex stiffness matrix which contains imaginary components representing 
the structural damping 

According to Reference 11, Chapter 6, the differential equation for this model is: 

T 8 2 u 9 2 u d 3u 

T = - u 6 — 

ax 2 at 2 at 


where 

T is the string tension (In this problem T is complex) 
y is the mass per unit length 
3 is the damping per unit length 

The finite difference representation for this equation is 



2u, 


Vr 


d 2 u. 


dt 
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The finite element model which corresponds to this equation is shown in Figure 2. Its 


equation is: 


m i u i + b i u i + (l+ig)k.(u i _ 1 - 2u i + u i+1 ) = P. 

where 

g = g 3 + g s 

is the structural damping. 

B. Input 


Parameters 


k i 

= 10 7 - 

scalar springs 

m i 

= 10.0 

scalar masses 

b i 

= 6.28318 - 

scalar dampers 

9 S 

= 0.05 

structural element damping 

^3 

= 0.05 

overall damping parameter 

N 

= 500 

number of scalar springs 


2. Constraints 

The end scalar springs are fixed on the outer ends so constraints are unnecessary. 

3. Eigenvalue Extraction Data 

Method : Determi nant 

Region of Interest: 5 < to < 16, -5.9 < a < 4.5 where p = a + iw ^ 

Normalization: Maximum deflection 


Answers: 

The natural frequency for an undamped string, according to Reference 11, is: 

t ml/T _ 7rnl/S 
w n £ ' N 


Its deflection shape is: 


or 


u(x) = 


sin 


mrx 


r in 


sin 


mn 

N 
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The modal masses are: 


l 0 0 m N 

M n = / UU Z (x)dx = 

0 


Substituting the real eigenvectors and eigenvalues into the complex equation for complex 
roots we obtain for each mode, n: 

H P 2 + (— p + (l+ig)a) 2 M = 0 
n 'm ' v y ' n n 


The solution is: 


P 



O+ig)^ 


A comparison of the complex roots is given in Table 
as the real eigenvectors, are nearly exact for the finite 


1. The eigenvectors, which are the same 
element model . 


Table 1. Comparison of NASTRAN and Analytic Complex Roots 


n 

Real Natural 
Frequency 

Theoretical Roots 
(radians per second) 

NASTRAN Roots 
(radians per second) 

■ 

1.0 

-.6285 ± 6.2753i 

-.6283 + 6.28321 

9 

2.0 

-.9425 ± 12.5621 

-.9419 + 12.5781 

9 

3.0 

-1.2566 ± 18.850i 
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RIGID FORMAT No. 7, Complex Eigenvalue Analysis - Direct Formulation 

Third Harmonic Complex Eigenvalue Analysis of a Gas-Filled Thin Elastic Cylinder (7-2-1) 
Fifth Harmonic Complex Eigenvalue Analysis of a Gas-Filled Thin Elastic Cylinder (7-2-2) 


A. Descri ption 

This problem demonstrates the use of symmetry to analyze specific harmonics of a fluid-filled 
structure. The problem to be solved consists of a cylindrical section filled with a compressible 
fluid. The end conditions for the cylinder and the fluid are two planes of antisymmetry , perpen- 
dicular to the axis. These end conditions correspond to the conditions that exist at periodic 
intervals along a long, fluid-filled pipe vibrating in one of its vibration modes. The antisym- 
metric boundary for the structure is defined by constraining the motions which lie in the plane. 

An antisymmetric boundary for the fluid corresponds to zero pressure. This may be modeled, in 
NASTRAN, by defining the plane of antisymmetry as a free surface with zero gravity. 

The lowest natural frequencies and mode shapes for the third and fifth harmonics are analyzed 
separately. For the third harmonic, the structure is defined by a section of a cylinder having an 
arc of 30 degrees or 1/12 of a circle. The fifth harmonic analysis uses a section having an arc 
of 18 degrees or 1/20 of a circle. The longitudinal edges, which were cut, are planes of symmetry 
and antisymmetry in order to model a quarter cosine wave length. 

The bulk data cards used are; AXIF, BDYLIST, CFLUID2 , CFLUID4 , C0RD2C, CQUAD1 , EIGC, FLSYM, 
FSLIST, GRIDS, MAT1 , PQUAD1 , RINGFL, and SPC1 . 


B. Input 

The finite element model for the third harmonic is shown in Figures 1 and 2. Parameters used 


are: 

8 = 2.88 x 10 3 lb/in 2 
pf = 1.8 x 10" 2 lb-sec 2 / in 4 

P s = 6.0 x 10' 2 lb-sec 2 /in 4 

E = 1.6 xlO 5 lb/in 2 

G = 6.0 x 10 4 lb/in 2 


(Bulk modulus of fluid) 

(Fluid mass density) 

(Structure mass density) 
(Young's modulus for structure) 
(Shear modulus for structure) 
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a 


10.0 inch 


(Radius of cylinder) 
(Length of cylinder) 
(Thickness of cylinder) 


£ = 10.0 inch 

h = 0.01 inch 

The model for the fifth harmonic is similar to the third harmonic model except that the angle 
covered by the structure is 18° instead of 30°. This is accomplished by simply removing the struc- 
tural elements and boundary GRIDB points corresponding to the two right-hand layers of structure 
(between 18° and 30°). The FLSYM, FSLIST and SPC1 cards are changed to reflect the 1/20 symmetry. 

C. Theory 

The derivations and results for this problem are described in Reference 16. The results for 
various dimensionless parameters are listed. The particular parameters for the problem at hand 
are: 

p f a 

n = p“FT 300 -° » ' 

/Gp f 

C = Vb^ = 2.5 . 

Pa 

ft = -gj- = 0.0 

where n is the ratio of fluid mass to structure mass. C is the ratio of the wave velocity in the 
structure material to the wave velocity in the fluid, ft is the factor describing static pressuri- 
zation, P . 

The basic assumptions for this analysis are: 

1. Thin shell theory is used for the structure. The bending moment terms in the force 
equilibrium equations are ignored in the results. 

2. The fluid is nonviscous, irrotational , and small motions are only considered. 

This particular problem becomes relatively easy to solve since the mode shapes for the fluid 
in a rigid container and the modes of the structure with no enclosed fluid have the same spatial 
function at the interface. Each mode of the fluid is excited by only one mode of the structure 
and each mode of the structure is excited by one mode of the fluid. The pressure in the fluid is 
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assumed to be a series of functions: 


p n e 1ut cos n<f> sin y-Q n (r,w) 


0) 


where Q n is a Bessel Function or a modified Bessel Function of the first kind. 
The characteristic shapes of the structure are a series of the form: 


u = A e la)t cos n<j> sin ~- 


( 2 ) 


where u is the displacement normal to the surface. The fundamental momentum equation for the 
fluid flow at the boundary is: 


V(P(r))*e 


p f u 


( 3 ) 


where e r is a unit vector in the radial direction. 

The forces on the structure at the boundary are: 


i g2 F i 

p (a) = -—2- - p s hu 


( 4 ) 


where the function is defined by the differential equation on the surface: 


w 4 F .Eh 3u 
' ' T ^ 


( 5 ) 


The solution for F-j is obtained by assuming that 


F-| = B e 1 *^ cos n(j) sin-— - 


( 6 ) 


Combining Equations 1 through 6 results in the relationships: 


2 „ _ p 3( V r ’“) 
P“ A - P n §7 


r=a 


( 7 ) 


Q n (a,co) P n 


a 2 Ah 


O IT2 + p s hw 


KM 


( 8 ) 
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Equation (7) is a statement of the continuity of displacement. Equation (8) states the 
ilance of the pressures. The above equations may be solved by iterating on oj. Reference 16 
rovides solutions for to over a wide range of parameters. 

0. Results 

The analytic and NASTRAN eigenvalues are listed in Table 1. The corresponding errors in the 
eigenvalues are tabulated and the maximum errors in displacement at the container wall are given 
as the percentage of the maximum value. The container displacements in the radial direction at 
4> = 0.0 are compared in Figure 3. 


Table 1. Comparison of analytical and NASTRAN results. 


Harmonic 

Natural 

Frequency (Hertz) 

Mode Shape 

Analytical 

NASTRAN 

% Error 

Max. % Error 
in Radial Displ. 

3 

1.579 

1.595 

1.0 

- 0.0 

5 

1.011 

1.049 

3.4 

0.5 
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Figure 3. Radial displacement for harmonic 3 and 5 normal modes. 



RIGID FORMAT No. 8, Frequency Response Analysis - Direct Formulation 

Frequency Response of a 10x10 Plate (8-1-1) 

Frequency Response of a 20x20 Plate (8-1-2) 

Frequency Response of a 10x10 Plate (INPUT, 8-1-3) 

Frequency Response of a 20x20 Plate (INPUT, 8-1-4) 

A. Description 

This problem illustrates the use of the direct method of determining structural response to 
steady-state sinusoidal loads. The applied load is given in terms of complex numbers which reflect 
the amplitudes and phases at each selected frequency. The steady-state response of the structure at 
each frequency is calculated in terms of complex numbers which reflect the magnitudes and phases of 
the results. Both configurations are duplicated via the INPUT module to generate the QUAD! elements. 

The particular model for this analysis is a square plate composed of quadrilateral plate 
elements as shown in 'Figure 1. The exterior edges are supported on hinged supports and symmetric 
boundaries are used along x = 0 and y = 0. The applied load is sinusoidally distributed over the 
panel and increases with respect to frequency. Although the applied load excites only the first 
mode, the direct formulation algorithm does not use this shortcut and solves the problem as though 
the load were completely general. 

B. Input 

1 . Parameters : 

a = b = 10 
t = 2.0 
E = 3.0 x 10 7 
v = 0.3 
y = 13.55715 

2. Loads: 

The frequency dependent pressure function is: 

P(x,y,f) = F(f) cos cos ^ 

where F(f) = 10- + 0.3f 


- length and width of quarter model 

- thickness 

- Young’s Modulus 

- Poisson's Ratio 

- nonstructural mass per area 
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3. 


Constraints : 


Only vertical motions and bending rotations are allowed. The exterior edges are hinged 
supports. The interior edges are planes of symmetry. This implies: 

along x = 0, 6^ = 0 

along y = 0, 6 x = 0 

along x = a, u z = 0 x = 0 

along y = b s u z = 8 y = 0 

all points, u x = u y " e z = 0 


D. Answers 

The excitation of the plate is orthogonal to the theoretical first mode. An explanation of 
the equations are given in Reference 8. The equations of response are: 


u z (f) 


F(f) 

(2ir) 2 p(f? - f 2 ) 


where f-| is the first natural frequency (10 cps). 

The following table gives the theoretical and NASTRAN results: 


Frequency 

u , x 10 4 
2 » 1 

cps 

Theory 

10x10 NASTRAN 

20x20 NASTRAN 

0 

1.868 

1 .874 


8 

6.435 

6.49 


9 

12.489 

12.69 

12.538 

10 

oo 

824.92 


11 

-11.833 

-11.67 
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RIGID FORMAT No. 9, Transient Analysis - Direct Formulation 
Transient Analysis with Direct Matrix Input (9-1-1) 


A. Description 

This problem demonstrates the capability of NASTRAN to perform transient analysis on a system 
having nonsyrrmetnc stiffness, damping and mass matrices. The problem also illustrates the use of 
time step changes, selection of printout intervals, application of loads, initial conditions, and 
a simple curve plot package. 

The matrices and loads used are actually the product of a transformation matrix and diagonal 
matrices. The resulting answers are easily calculated while the input matrices are of general 
form. The matrix equation solved is: 

[M]{ If) + [B]{ u} + [K] { u } = {P(t)J 


The problem is actually four disjoint single degree of freedom problems which have been 
transformed to a general matrix problem. Figure 1 illustrates the problems schematically. 


The resulting diagonal matrices are premultiplied by the matrix: 


[X] = 


2 

-1 

0 

0 


-10 0 
2-10 
-1 2 -1 
0-12 


The answers for the disjoint problem above will be the same as for the general matrix problem 
since the general case: 


[X]([M 0 ]{UJ + CB 0 ]{u> + CK 0 ]){u} = [X]{P} 


has the same results as the disjoint case: 


[M 0 ]{uJ + [B q ] { u} + [K q ] { u } = {P} 
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B. Input 



The 

initial 

condi ti ons 

are: 


u 10 

= 0 


fl 10 

= 10.0 

U 11 

= 0 


“n 

= 0.5 

u 12 

= 0 

• 

u 12 

= 0 

u 13 

= -10, 

.0 

<13 

= 0 


3. At t = 1.0 a step load is applied to each point. The load on the uncoupled problems is: 
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The transformed load is: 


{P} = [X]{P 0 > 


-1.5 

- 1.0 

-13.5 

36.0 


C. Answers 

The results are responses of single degree of freedom systems. Equations are given in 
Reference 12, Chapter 9. 

0 < t < 1.0 , At = .005 


U 10 

= sin lOt 



u 10 

U 11 

= 0.05(1 - 

e‘ 10t ) 


& 11 

u 12 

= 0 



Q 12 

U 13 

= -10e -10t 



*13 

t > 1.0 , 

At = .015 



u 10 

= sin lot 




U 11 

= 0.05(1 - 

e-10t, 

+ 0.1 (t ■ 

- 1.1 

u 12 

= 0.04 (l - 

e” 3t [cos4(t-l ) + 

|si 

u 13 

■ -10e- ,ot 

+ 1 


0 


-10t 


-lot 


-10(t-l) 


Figures 2 through 5 are tracings of the NASTRAN plots of the functions. The deviations of 
the NASTRAN results and the theoretical response are due to the selection of time steps. For 
instance point 11 has a time constant equal to two time steps. The initial error in velocity due 
to the first step causes the displacement error to accumulate. Using a smaller time step has 
resulted in much better results. 
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9 . 1-4 






Displacement, u 


NASTRAN 

O Calculated 



Figure 2. Point 10, displacement. 
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Displacement, u 



Figure 3. Point 11, displacement. 
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0.036 
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0.020 

0.016 
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0.008 

0.004 

0 






Displacement, u 



Figure 5. Point 13, displacement. 
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RIGID FORMAT No. 9, Transient Analysis - Direct Formulation 

Transient Analyisis of a 1000 Cell String, Traveling Wave Problem (9-2-1) 

Transient Analyses of a 1000 Cell String, Traveling Wave Problem (INPUT, 9-2-2) 

A. Descn ption 

This problem illustrates the ability of NASTRAN to perform time integration studies using 
the structural matrices directly. At each time step the applied loads, the structural matrices, 
and the previous displacements are used to calculate a new set of displacements, velocities, and 
accelerations. Initial displacements and velocities are also allowed for all unconstrained coor- 
dinates. The INPUT module is used to generate the scalar springs and masses. 

The structural model consists of a 1000 cell string under constant tension modeled by scalar 
elements. The string is given an initial condition at one end consisting of a triangular shaped 
set of initial displacements. The wave will then travel along the string-; retaining its initial 
shape. The ends of the string are fixed causing the wave to reflect with a sign reversal. 

Figure 1 illustrates the problem and the scalar element model for each finite increment of 
length. 

B. Input 

1 . Parameters : 

k i = Ax' = ~ scalar s P nn 9 rates 

m^ = pAx = 10 - scalar masses 

N = 1000 - number of cells 

where 

T is the tension 
Ax is the incremental length 
p is the mass per unit length 
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2. Loads: 


The initial displacements are; 


u 2 = 

.2 

u 12 

= 1.8 

u 3 = 

.4 

u 13 

= 1.6 

u 4 = 

.6 


* 

'll = 

2.0 

U 21 

= 0.0 


C. Answers 


As shown in Reference 11. Chapter 6, the wave velocity c is, 


c = 




±1000 points/umt time 


The imti 
The first wave 
a sign change, 
traveled their 


al displacement may be divided into two waves, traveling in opposite directions, 
travels outward; the second wave travels toward the fixed support and reflects with 
The theoretical and NASTRAN results are compared in Figure 2, when both waves have 
complete width. 
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3 


i- 


\ 


T 


> 

P = Hass/Length 


u 999 u 1000 

iX 


1000 Cell String 


1-1 


J i+1 


-o O W\ — o- 

m i-l V m i k i+l m 1+l 


Finite Element Model 


Figure 1. Representations of dynamic string. 













RIGID FORMAT No. 9, Transient Analysis - Direct Formulation 
Transient Analysis of a Fluid-Filled Elastic Cylinder (9-3-1) 


A. Description 

The fluid-filled shell, used for analysis of the third harmonic, in Demonstration Problem 
No. 7-2 is subjected to a step change in external pressure at t = 0 of the form 

. 7TZ 

p = P 0 Sin — cos n<f> 

The fluid is assumed incompressible in order to obtain an analytical solution with reasonable 
effort. The harmonic used is n = 3. 

In addition to the cards of Demonstration Problem No. 7-2, DAREA, PRESPT, TL0AD2, and TSTEP 
cards are also used. Selected displacements and pressures are plotted against time. 


B. Input 

The finite element model is shown 

B = oo 

p f = 1.8 x lCf 2 lb-sec 2 /in 4 
P s = 6.0 x 10“ 2 lb-sec 2 /in 4 

E = 1.6 x 10 5 lb/in 2 

G = 6.0 x 10 4 lb/in 2 
a = 10.0 inch 

£ = 10.0 inch 
h = 0.01 inch 



n Figures 1 and 2. Parameters used are: 

(Bulk modulus of fluid - incompressible) 
(Fluid mass density) 

(Structure mass density) 

(Young's modulus for structure) 

(Shear modulus for structure) 

(Radius of cylinder) 

(Length of cylinder) 

(Thickness of cylinder wall) 

(Pressure load coefficient) 
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C. Theory 


The theory was derived with the aid of Reference 16 as in Demonstration Problem No. 7-2. 
Since the fluid is incompressible, it acts on the structure like a pure mass. Neglecting the 
bending stiffness, the equation of force on the structure is: 


1 

P- = (m + m f ) w + -~ 
s T a dz d 


( 1 ) 


where: 

P s is the loading pressure on the structure (positive outward), 
m - p s h is the mass per area of the structure, 
ny is the apparent mass of the fluid, 
w is the normal displacement (positive outward) 


The function F is defined by the equation, 


v 4 f 


Eh 

a 


a 2 w 


(2) 


The spatial functions of pressure, displacement, and function F may be written in the form 


P$ 55 P 0 sin — cos n<J> 


TTZ 

w = w 0 sin cos n <j> 


r r- . TTZ , 

F = F sin -y- cos n4> 

0 As 


where p 0? w Q > and F Q are variables with respect to time only. 
Substituting Equations 3 into Equation 2 we obtain: 


w ° 


[> * (r 12 


jy 


Substituting Equations 3 and 4 into Equation 1 we obtain: 


(3) 
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(m + m f ) w Q + 


Eh 


(5; 


* e) i 


2T2 w o 


The incompressible fluid is described by the differential equation: 


V p = 0 


( 6 ) 


Applying the appropriate boundary conditions to Equation 6 results in the pressure distribu- 


tion: 


p = p r sin f cos (nt}>) I p (2f) 


(7) 


where I n is the modified Bessel function of the first kind and p r is an undetermined variable. The 
balance of pressure and flow at the boundary of the fluid, with no structural effects, is described 
by the equations: 


Pr Uf) 


( 8 ) 


p f" = 'If- 


r=a 


(9) 


Substituting Equations 3 and 7 into Equation 9 results in: 


- tt t 1 ( ffa \ n 

Pf w o ' n nu' P r 


( 10 ) 


Eliminating P with Equations 8 and 10 gives the expression for apparent mass, m^: 


_ T /ira\ p f w o 

p ° ‘ " (t) ^ 
l n\ l / 


= m f w o 


01 ) 


Substituting the expression for ny from Equation 11 into Equation 5 results in a simple single 
degree of freedom system. When the applied loading pressure is a step function at t = 0, 
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where 


w 


(1 - cos tot) sin— cos n4> 


02 ) 


!!a 

k 


anG 


and 



k 


Eh 

a2 [’ * ©T 


I (™) 

A U A SL n\£ / 

rrv = m + rrv = p h + — — r— ■ 

T f S T 7T T » /TTa \ 

n\ t ) 


D. Resul ts 

A transient analysis was performed for the case n = 3 on the model and various displacements 
and pressures were output versus time up to one second. The theoretical frequency is calculated to 
be 1.530 Hertz and the period is 0.633 seconds. The displacements at two points on the structure 
(Point 91 is located at <j> = 0, z = 5.0; Point 94 is located at <f> = 18°, z = 5.0) are plotted 
versus time in Figure 3. 

The maximum error for the first full cycle occurs at the end of the cycle. The ratio of the 
error to maximum displacement is 4.75%. Changes in the time step used in the transient integration 
algorithm did not affect the accuracy to any great extent. The most probable causes for error were 
the mesh size of the model and the method used to apply the distributed load. The applied load was 
calculated by multiplying the pressure value at the point by an associated area. The "consistent 
method" of assuming a cubic polynomial displacement and integrating would eliminate the extraneous 
response of higher modes. The method chosen in this problem, however, is typical of actual 
applications. 
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cure 2. Transient analysis model. 



Theory 


o 
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I 


i 
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Figure 3. Displacement at midpoint (z = 5.0). 



RIGID FORMAT No. 9 (APP HEAT), Linear Transient Heat Transfer Analysis 
Plate with Suddenly Applied Flux and Edge Temperature (9-4-1) 


A. Description 

The time history of the temperature in a lonq thin plate initially at zero deqrees is analyzed 
using NASTRAN's transient heat analysis capability. At time t=0 a heat flux is applied on one sur- 
face of the plate and simultaneously the temperature along the edges is increased. These tempera- 
tures are maintained at a value by using a large heat flux through a good conductor to ground. The 
problem is one dimensional since it is assumed that no temperature variation exists along the length 
or through the thickness. Since the plate is symmetric about the center plane, only one half of the 
plate is modeled. 

B. Input 

The plate is shown in Figure 1 and the idealized NASTRAN model, shown in Figure 2, is repre- 
sented by five R0D elements going from the centerplane to the edge. The conductor-ground arrange- 
ment is modeled by an ELAS2 element and an SPC card referenced in Case Control. The injected heat 
flux at the edge is specified using DAREA and TL0AD2 cards which are referenced in Case Control 
through a DL0AD card. The surface heat flux is specified on a QBDY1 card and references the TL0AD2 
card. The time step intervals at which the solution is generated are given on the TSTEP card. The 
initial temperature conditions are specified on the TEMPO card and referenced in Case Control by an 
IC card. The heat capacity and conductivity are given on the MAT4 card. 

C. Theory 

The analytic solution is 

0.5 |^1 - J e ' (2n+1)2t cos(2n+l W2j + 

50. [(,-„■) - ^ £ (2n+T)~3 e "^ 2n+1 ^ ^os^n+l )nx/2j 

D. Results 

A comparison of theoretical and NASTRAN results is given in Table 1. 
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Table 1. Theoretical and NASTRAN temperatures. 



GRID(X) 


10(0.) 

12(.2) 

1 4 ( . 4 ) 

1 6( . 6 ) 

18(.8) 

20(1.) 

Theory* 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

t = 0 







NASTRAN 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

Theory* 

31.282 

30.222 

26.952 

21.204 

12.562 

.500 

t = 1 







NASTRAN 

30.641 

29.612 

26.433 

20.826 

12.362 

.500 

Theory* 

43.430 

41.776 

36.780 

28.344 

16.316 

.500 

t = 2 

i 






NASTRAN 

43.117 

41.478 

36.527 

28.160 

16.218 

.500 

Theory* 

47.916 

46.026 

40.396 

30.971 

17.696 

.500 

t = 3 







NASTRAN 

47.755 

45.890 

40.280 

30.887 

17.652 

.500 

v * °o Theory 

50.500 

48.500 

42.500 

32.500 

18.500 

.500 


* 


n = 0 term only. 
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Figure 2. Idealized NASTRAN model. 
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RIGID FORMAT No. TO, Complex Eigenvalue Analysis - Modal Formulation 
Rocket Guidance and Control Problem (10-1-1) 


A. Description 

This problem, although a simplified model, contains all of the elements used in a linear 
control system analysis. The flexible structure, shown in Figure 1, consists of three sections: 
two sections are constructed of structural finite elements; the third section is formulated in 
terms of its modal coordinates. A sensor is located at an arbitrary point on the structure and 
connected to a structural point with multipoint constraints. The measured attitude and position 
of the sensor point is used to generate a control voltage for the gimbal angle of the thrust 
nozzle. The nozzle control is in itself a servomechanism consisting of an amplifier, a motor, 
and a position and velocity feedback control. The nozzle produces a force on the structure due 
to its mass and the angle of thrust. The motion of any point on the structure is dependent on 
the elastic motions, free-body motions, and large angle effects due to free-body rotation. 

The guidance and control system is shown in block diagram form in Figure 2. The definitions 
for the variables and coefficients along with values for the coefficients are given in Table 1. 
The use of the Transfer Function data card (TF) allows the direct definition of the various rela- 
tions as shown in Figure 2. 

i 

B. Modeling Techniques 


1. A section of the structure is defined by its modal coordinates by using a modification 
of the method given in the NASTRAN Theoretical Manual. The algorithm is given as 


follows: ** 

Define £. , i = 1 , n - modal deflections scalar points 

u r - grid point components used as nonredundant supports for modal 
test. These may or may not be connected to the rest of the 
structure . 


u c - grid point components to be connected to the remaining 
structure (not u r points) 

x , i = 1, n - rigid body component degrees of freedom for the nonzero modes 
The relations between these variables are defined by using.mul tipoint constrai nts «wi/th 


the following relationships* 
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a) 


l"c> ■ [“cr^V 

<|> cl . is the angular deflection of point u c for mode i. D cr is the deflection of 
point u c when the structure is rigid and point u r is given a unit deflection. 

b) {xp = DC 1 3“ 1 [H] T {u r > = [G]{u r ) 

[K.] is a diagonal matrix. Each term K , the modal stiffness, is defined as: 

K-, = f o) 

where nu is the modal mass and w is the natural frequency in radians per second. [H] 
is determined by the forces on the support points due to each nonzero eigenvector 

P r ■ -l H n 5, <«, * 0) 

c) Scalar masses and springs are connected to each modal coordinate as shown by 
Figure 3a. 

d) The structure to be added in this problem consists of a simply supported uniform 

beam as shown in Figure 3b. The support points, u^, are y^ and y^g. The additional 
degree of freedom to be connected is u c = 0^. Four modes are used in the test 

problem. The following data is used to define and connect the modal coordinates of this 
substructure. 

The mode shapes are: 


mrx 


♦n (x) ' sin Z 


The modal frequencies, masses, and stiffness in terms of normal beam terminology are: 

n = 1 , 2, 3, 4 


2 2 

n V El 


pA 


pA& 

2 


4 4 

n tt El 


2r 


The forces of support for each mode are: 
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P y 06) = I - " 

y SL 

PJ19) = [H) n ^ n 3 

y l 

The motion is defined by multipoint constraints: 

0 16 = I ( y 19 ' y 16 ^ + l IT 

n 

The free-body components of the modes are defined, using multipoint constraints, as: 



2. The mass of the nozzle would normally be included with the structural modeling. How- 
ever, to demonstrate the flexibility of the Transfer Function data, it is modeled as part of 
the guidance system as shown in Figure 4a. 

Defining the angle of thrust, y, to be measured relative to the deformed structure, the 
forces which result are: 

T = (I no + x n m i )( ? + * e V “ m n X n y l 

F y = m n y l - x n m n (y + h \ > " F n y 
Using the thrust force, F , as a constant, the transfer functions are: 

! n y * T + r n $2 9 1 - x n m n ^ y l = 0 

m n $2 y l - (x n m n ^ + F n )y ' x n m n s * 9 1 = 0 

(0)e 1 + t = o 
2 

where: I = I ^ x m = 500 

n no n n 
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3. The large angle motion must be included in the analysis since it contributes to the 
linear terms. The equations of motion of the structure are formed relative to a coordinate 
system parallel to the body. As shown in Figure ^b, the accelerations are coupled when the 
body rotates . 

Since the axial acceleration, x, is constant throughout the body, the vertical acceler- 
ation at any point, to the first order, is: 

y abs = y rel + x6 l = y rel + y e 

An extra degree of freedom y 0 is added to the problem and coupled by the equations: 

my e = F n e ] 

y abs ~ y rel + y 9 

4. The center of gravity (point 101) and the sensor location (point 100) are rigidly con- 
nected to the nearest structural point with multipoint constraints. For instance the 
sensor point is located a distance of 4.91 from point 8 as shown in Figure 4c. 

It is desired to leave point 101 as an independent variable point. Therefore point 8 
is defined in terms of point 101 by the equations: 

y 8 = y 101 + 4 ‘ 91e l01 

0 8 = 6 101 


C. Answers 

A comparison of the NASTRAN complex roots and those derived by a conventional analysis 
described below are given in Table 2. The resulting eigenvectors were substituted into the 
equations of motion to check their validity. The equations of motion for a polynomial solution 
may be written in terms of the rigid body motions of the center of gravity plus the modal dis- 
placements. The equations of motion using Laplace transforms are: 

m$2 y cg = F n (6 l + y) 

I$2 6 eg = ” F n X 1 Y 
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The inertia forces of the nozzle on the structure may be ignored. 
The motion of the nozzle, as explained in section B-2, is: 


(j- + ts + l^y = (a + bs)y $ + (c + ds)e s - j- 0, + 


s^m x 


n n 


yi 


( y is defined as the relative angle between the nozzle and the structure.) 

The flexible motions at the sensor point, y $ and e s , may be defined in terms of the modal 
coefficients and the rigid motions of the center of gravity. 


*s 

6 

S 


0 

eg 


+ * 2 e cg + l 00 , l ? i 


+ 


l *. 


100,i 


5 

1 


The motions of the nozzle point, in terms of the modal and center of gravity motions are: 


0 


1 


0 + 
eg 



S 

i 


y, = y X-. 0 + I <l> C 

J 1 J cg 1 eg ^ l,i i 

The modal displacements are due primarily to the vertical component of the nozzle force. Their 
equation of motion is: 

= F n v 


where 

6. * is the deflection of point j for mode i 
J 

i 

4>- . is the rotation of point j for mode 1 
J j i 

m 1 is the modal mass of mode i 

is the natural frequency of mode i 
is the modal displacement of mode i 

Using two flexible modes the characteristic matrix of the problem is given in Figure 5. The 
determinant of the matrix forms a polynomial of order 10. The roots of this polynomial were 
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,sd by a standard computer library routine and are presented in Table 2 as the analytical 
Its. The rigid body solution is also presented. 

The differences between the two sets of answers is due to the differences in models. The 
TRAN model produces errors due to the finite difference approximation and the number of modes 
oen to model the third stage. The polynomial solution produces errors due to the approxima- 
tes used in the,equations of motion as applied to control system problems. 

As a further check the first eigenvalue (A = -7.41) was substituted intoj the matrix given m 
Figure 5 and the matrix was normalized by dividing each row by its diagonal value. The NASTRAN 
eigenvector was multiplied by the matrix, resulting in an error vector which theoretically should 
be zero. Dividing each term in the error vector by its corresponding term in the eigenvector 
resulted in very small error ratios. 
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Table 1. Variables and Parameters 


Extra Point 


Nianber 

Symbol 

Description 

1010 

e y 

Voltage describing y 

- 1011 

e e 

Voltage describing e 

1020 

E yc 

Control voltage for y (Input) 

1021 

E 0C 

Control voltage for q (Input) 

1030 

E 

T 

Attitude error function 

1040 

£ 

y 

Nozzle position error 

1050 

E m 

Voltage for Nozzle servo 

1060 

T 

Torque for Nozzle servo 

1070 

y 

Nozzle Thrust angle relative to structure 

1080 

h 

Position increment due to attitude 


Parameters 

Value 

Descri ption 

K s 

1.0 

Servo ampl l fi er gai n 

•V. 

500 

Servo gain 

T 

.1414 

Nozzle angular velocity feedback 

x n 

3.0 

Distance from nozzle C.G. to Gimbal axis 

Z n 

500.0 

Inertia of Nozzle about gimbal axis 

F n 

4.25x1 0 6 

Thrust of Nozzle 

m n 

50 

Nozzle mass 

B e 

100.0 

Overall vol tage-to-angle ratio 

6 y 

1.0 

Overall voltage to postton ratio 

a 

.16 

Position feedback coefficient 

b 

.28 

Velocity feedback coefficient 

c 

15.0 

Angle feedback coefficient 

d 

7.0 

Angular velocity feedback coefficient 

m 

8.5xl0 4 

Mass of structure ^ 
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Table 2. Comparison of Complex Roots for NASTRAN Modeling vs. Simplified Polynomial Expansion 


Rigid Body Model 

2 Flexible Modes Model 



NASTRAN* 

POLYNOMIAL 

NASTRAN 

POLYNOMIAL 

- 540 ± .8211 

-.522 ± .8021 

-.507 ± .8191 

-.494 ± .8011 

-1.68 ± Oi 

-1.74 ± Oi 

-1.41 ± Oi 

-1 .46 ± Oi 

+.751 ± 5.96i 

+.774 ± 5.98i 

+.520 ± 3.82i 

+.522 ± 3.83i 


*Not published. 
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Structure 



a) Diagram for input of modal data 



Figure 3. Modal data input diagrams. 
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RIGID FORMAT No. 10 (APP AER0) , Aeroelastic Analysis 
Aeroelastic Flutter Analysis of a 15° Swept Wing (10-2-1) 


A. Description 

This problem illustrates the use of the aeroelastic analysis to determine flutter frequencies 
and mode shapes for an untapered wing having 15° sweep and an aspect ratio of 5.34 as shown in 
Figure 1. 

B. Input 

Bulk data cards used include CAER01 , PAER01 , SPLINE2, SET1, AER0, MKAER01 , FLUTTER, and FLFACT 
as illustrated in User's Manual Section 1.11. 


C. Theory 

Reference 22 specifies the reduced frequency k = .1314 (p.17), frequency ratio <o/u) a = 0.51 
(p.35) and torsion frequency to a = 1488 (p.17). 

The flutter velocity is found from 

REFC x w x to 

— O — a — 

.2 a) 

v = y - = £ = 5980 m/sec. 


where REFC is the reference length input on the AER0 bulk data card. 
The flutter frequency is found from 


f = 


0) 


a 



2tt 


= 121 Hz 


D. Results 

The results obtained are compared with both theoretical results using the modified strip 
analysis method and with experimental results. The flutter velocity is in good agreement. (See 
Figure 2. ) 

Frequencies are automatically output while mode shapes used in the modal formulation are 
obtained using an ALTER to the Rigid Format following the Real Eigenvalue Analysis Module. 

Mode shapes for all points in the model may be obtained by checkpointing the problem using the 
Normal Mode Analysis (Rigid Format 3) and subsequently restarting using the Aeroelastic Analysis. 
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RIGID FORMAT No. 11, Frequency Response - Modal Analysis 
Frequency Response and Random Analysis of a Ten-Cell Beam (11-1-1) 

A. Descri ption 

This problem demonstrates the frequency response solution of a structure using uncoupled 
modal formulation. With modal formulation, the structural degrees of freedom used in the solu- 
tion are the uncoupled modal displacements. The solution equations are simple and efficient. 

The saving in time, however, is offset by the operations necessary to extract the modes, transform 
the loads to modal coordinates, and transform the modal displacements to structural displacements. 

This problem also illustrates the various methods of applying frequency response loads. 

Loads may be input as complex numbers, with phase lag angles and/or time lag factors. The loads 
may be added together for each subcase. 

The structure to be solved consists of a beam with simple supports on the end as described 
in Figure 1. The parameters selected produce natural frequencies of 50, 200, 450 and 800 cps. 

The applied loads for the three subcases are applied to the center with variations in phase angles, 
time lags and input formats. The first two subcases use three loaded points which, in essence, 
simulate a load on the center. 

Included in the structural representation is a "general element" representing the first two 
cells of the ten-cell beam. The flexibility matrix, [Z], of the element represents the displace- 
ments of grid points 2 and 3 when point 1 is fixed. The rigid body matrix, [S], represents the 
rigid body motions of points 2 and 3 when point 1 is displaced in the x, z, or directions. 

The random analysis data consists of a flat power spectral density function ("white noise") 
for the three loading subcases. The first subcase spectral density is connected to the third 
subcase spectral density, simulating two interdependent probability functions. The XY-plotter is 
used to plot the displacement and acceleration power spectral density function of grid 6 (center 
of the beam). The displacement autocorrelation function is also plotted for the same point. All 
values are tabulated in the printout. The NASTRAN power spectural density results are compared 
against a simplified analytic calculation in Figure 2. 

A static analysis restart of the frequency response problem is demonstrated. Gravity and 
element enforced deformation loads are used with a change in the single-point constraints. 
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B. Input 


1 . Parameters : 

1 = 20 
^ = .083 

A = 21.18922 

E = 10.4 x 10 6 

p = .2523 x 10 

M = pAX, 


- length 

- bending inertia 

- cross sectional area 

- modulus of elasticity 

3 

- mass density 

- total mass 


2. Constraints: 

Uy = e x = e z = ° 

u xl = u zl = u zl 1 = 0 

u xl = u zl = u xll = u zll 


- all points 

- frequency response 
0 - static analysis 


3. Modal Data: 

Interval: 40 < f < 1000 cps 

Normalization: Modal Mass = 1.0 

Number of modes used in formulation: 4 

-4 

Modal Damping ratio: g = 4 x 10 f 

4. Loads, Frequency Response: 

The loading functions for subcase 1 are: 


P z,5 

= 50 


M y,5 

= -100 


P Z,6 

= 50 + 

100(cos60° + i sin60°) 

^ set'IT^ ^ 

P z,7 

= 50 


M y>7 

= 100 
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The loading for subcase 2 is: 


P z,5 

50 



M y,5 = 

-100 



P z,6 = 

50 + 

ig0(cos2f° ^ 

l sin2f°) 



SET 7, t = 

.005555 

P z,7 " 

50 



M y,7 = 

100 




The load for subcase 3 is: 

P 7 , = 2[75 + 50i(cos30° - i sin30°)] = 200 + 86. 6i 
z ,o 

Note: At f = 30cps the three subcases are nearly identical. 

5. Random Analysis Data 

The nonzero factors for the three subcases are: 


$ u = 50 


H3 


S 31 - 50 


0 < f < 100 


S 22 = 100 
S 33 = 50 


s = 0 , f > 100 

* J 

The time lags selected for the autocorrelation function calculations are: 


t = 0 . 0 , 0.001 , 0 . 002 , 0.1 


6. Static Loads for Restart 

The problem is run first as a frequency response analysis. It is restarted as a static 
analysis with the following loads: 

Gravity vector: g z = 32.2 

Element Deformation: = 0.089045 (expansion) 
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C. Analysis 


1. The theoretical eigenvalue data, according to Reference 8 is: 


(2ttH 


El 

2 V A 


T- = 50, 200, 450, 800 


(natural frequencies) 


m n = 1.0 (modal mass) 


♦„(X) = [/ pA sin 2 ^ dxj 1 sin(^) sin ( I X^) ( mode shape) 


1 


2. The theoretical frequency response at the center point is essentially the response of 
the first mode which is: 

\ *1,6 P j (u) *1J 

Ug(w) = — * (j = degree of freedom number) 

m-j (o)^ - u) 2 + igaxi)^ ) 

At the first natural frequency of 50 cps, the response will be nearly equal to the 
response of the first mode. The response at the center point for the three subcases 
are: 


Subcase 1 and 3 


..1 3 94.764 + 41. 033i 

u = u = 

° 6 (50-f 2 ) + if 


The results are: 


f 

0 

30 

50 


Ug (one mode) 

.0413 @ 23.42° 
.0646 @ 22.34° 
2.066 0 293.42° 


u 6 (NASTRAN) 

.0429 0 22.9° 
.0668 @ 21.8° 
2.074 0 281.5° 
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Subcase 2 


u 2 = 23.691 (3 -i- 2cos2f - 2i sin2f) 

6 (50 - f 2 ) + if 

Theoretical and NASTRAN results are: 


f u & (one mode) u g (NASTRAN) 


0 

30 

50 


.047 @ 0° 

.0646 @ -22.34° 
1.565 @233.4° 


.049 @ 0° 

.0668 @ -23.97° 
1.577 @ 223.0° 


3. The random analysis is explained in Reference 15. The power spectral response coeffic- 
ients for the three subcases are given by the matrix: 




0.5 



0 0.5 

1.0 0 

0 0.5 


If (H } is the vector of the responses of a point, j, to the three loading cases, the 

J 

power spectral response, S , is: 

S = {H -) T [S ]{H } (H. is the complex conjugate) 

J J J J 

or 

Sj = lOOCO.SlH^.I 2 + H 3j + H 3j H^) + | H 2 J 1 2 + 0.5|H 3j | 2 ] 


Since H^. 


j , then : 


Sj = 200|H lj | 2 + 100|H 2j .| 2 


The mean square response in obtained by integrating the power spectral density over the 
frequency. In this particular case the frequency increments are uniform and the mean 
square response is simply 


E 


i 



- Sjtf^lAf 
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The analytic solution for the displacement spectral ^ dens i ty response of the center point 
due to the first mode is: 


s = 200(1. 066xl0 4 ) + 100(.5613x10 3 )(13 + 12cos2f) = 2.862x10 6 + ,6735x10 6 cos2f 
J [(50 2 - f 2 ) 2 + f 2 ] (f 4 - 4999f 2 + 50 4 ) 

The results of the above equation are compared with the NASTRAN results in Figure 2. 

The mean deviation, a , is: 

J 

-J r i 

°d - f 0 ) 

where f p and f Q are the upper and lower frequency limits, was checked by summing the - 
NASTRAN results. 


4. The results of the static analysis restart are: 

a) The gravity load produces normal displacements (in the z direction) and element 
moments as follows: 

u z ( x) = fjjjf (l* - 2 Ax 2 + x 3 ) 


M ] (x) = ^2- (x 2 - Ax) 


b) The element deformation produces the following axial forces and displacements: 

u x = ■ AT x (x < 18) 

In numerical terms the displacements of the center point (x = 4) are: 

Theoretical NASTRAN 


'x 6 

= 4.452 x 

C\J 

0 

4.435 

x 10" 2 

26 

= 4.155 x 

0 

1 

4.121 

x 10' 4 
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The element forces at the center of the beam are: 


Theoretical 
F x5 = -.9811 x 10 6 


NASTRAN 
-.9848 x 10 6 


M 6 = -8.607 


-8.607 
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S 2 (f), Point 6 



Figure 2. Power spectral density of center point displacement. 
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RIGID FORMAT No. 11, Frequency Response Analysis - Modal Formulation 

Frequency Response of a 500-Cell String (11-2-1) 

Frequency Response of a 500-Cell String (INPUT, 11-2-2) 


A. Description 

This problem illustrates the solution of a large frequency response problem using modal coordi- 
nates. When large numbers of frequency steps are used, or the problem is very large, the relative 
efficiency of the modal formulation is more attractive than the direct formulation. The structural 
model consists of scalar points, springs, and masses which simulate the transverse motions of a 
string under tension, T, with a mass per length of y. The model and its finite element representa- 
tion is shown in Figure 1. A duplicate model is obtained via the INPUT module to generate the scalar 
springs and masses. 

Selected scalar point displacements and scalar element forces are plotted versus frequency. 

The magnitude and phase of the displacements are plotted separately, each on one-half of the 
plotter frame. The magnitude plots for the selected points are all drawn on a whole plotter frame 
for comparisons. The center spring element has the magnitude of its internal force plotted versu c 
frequency. 

B. Input 

1. Parameters* * ‘ 

m. = 10 - mass 

K. = 10 7 - spring rate 

N = 500 - number of cells 

where 

K i = ST * m i = ^ Ax 

2. Loads 

The load on each point is: 

P^w) = Axp x = IOtt 3 

where p x is the load per length of string. 
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The steady state frequency response is desired from .1 to 10 cycles per second in 15 
logri thmiC) increments. 


3. Real Eigenvalue Data 

Method: Inverse Power 

Region of interest: 0 < f < 21 

Normalization: Mass 

Number of modes used in formulation: 20 


C. Answe rs 

The analysis of the string is given in Reference 11, Chapter 6. The response, £ , of mode 
number n is given by the equation: 


/ P(x) sin(^) dx 


J ~J l psin 2 (~) dx J [u> 2 - w 2 ] 


/K 

where co n , the natural frequencies, are for the theoretical continuous string. 


For a uniform Load: 


/ Hx) Sln^~) dx 


mr 


2P.N 


10V 


/psin 2 (H*) dx 


2 


Nm i 


2.5 x 10° 


The displacement of the center point is: 


u(f) = I C n sin f- 




For instance at f = 0.1 the response due to 20 modes is: 

u(|) = .97895 (Theory) 

u 251 = .97888 (NASTRAN) 
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RIGID FORMAT No. 12, Transient Analysis - Modal Formulation 
Transient Analysis of a Free One Hundred Cell Beam (12-1-1) 


A. Description 

The problem demonstrates the transient analysis of a free-body using the integration algorithm 
for uncoupled modal formulations. The model is a hundred-cell beam with a very large mass attached 
to one end as shown in Figure 1. Modal damping is included as a function of natural frequency. It 
does not affect the free-body (zero frequency) modes. The omitted coordinate feature was used to 
reduce the analysis set of displacements to correspond to eleven grid points. 

Both structure plots and curve plots are requested. The types are as follows: 

1. Stereoscopic structure plots of the deformed structure are drawn for a specified time 

step. 

2. Orthographic projections of the deformed structure are plotted. However, two variations 

are plotted on each frame. The bottom region of the frame shows the deformed shape and 

the top region shows vectors at every tenth grid point which are proportional to the z- 

displacement at each specified time step. 

3. Curve plots and printout of displacement versus time and of acceleration versus time 
are requested. 

When a structure is used without additional transfer functions or direct matrix inputs,' the 
transient analysis solves exact equations for the uncoupled modes. . The only errors will be in 
the discarded modes and the straight line approximation of the loads between time ‘Steps; The 

speed of this solution is offset by the fact that the eigenvalue calculation is relatively costly 

» 1 > 

and the transformation of the vectors to and from modal coordinates could be time consuming. 

The mass and inertia on point (1) were selected to be much larger than values of the beam. 

The answers will therefore approximate a beam with a fixed end. 
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B. 


Input 


1 . Parameters 


Beam: 

% 

I 

A 

E 

P 


20 

.083 

1.0 

10.4 x 10 6 
.2523 x 10' 3 


(Length) 

(Bending inertia) 

(Cross sectional area) 
(Modulus of elasticity) 
(Mass density) 


Lumped Mass: 


m^ = 10.0, I 22 ] = 1666.66 


2. Damping: 

The damping coefficient for each mode is a function of the natural frequency. The 
function is: 


g = 10‘ 3 f 


3. Load: 

P ir ._ = 100 sin(2ir*60t) 

z , 1 01 

4. Real Eigenvalue Data 

Method: Inverse Power 

Region of Interest:' 0 < f < 1000 
Normalization: Mass 

D. Answers 

The NASTRAN results are compared in Figure 2 to the analytic results which use one mode. The 
modal mass may be calculated using the formula for the mode shape given in Reference 8. The modal 
displacement is a single degree of freedom response with a closed form solution. 
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RIGID FORMAT No. 13, Normal Modes with Differential Stiffness 
Normal Modes of a 100-Cell Beam with Differential Stiffness (13-1-1) 


A. Description 

This problem illustrates the effects of differential stiffness on the solution for the normal 
modes of a beam under axial compression. 

The natural frequencies of the beam are affected by this load as shown in Reference 23. The 
loading specified here is one half of the Euler value for compression buckling which decreases the 
unloaded natural frequency, w, proportional to 


7T 2 El 
Z 2 



where F is the applied load. 

The structural model illustrated in Figure 1 is a uniform 100 cell beam hinged at both ends. 


B. Input 

1. Parameters: 


A = 2.0 
I = 0.667 
E = 10 4xl0 6 
Z = 100.0 
p = 2. 0x10" 4 


(cross sectional area) 
(bending inertia) 
(modulus of elasticity) 
(length) 

(mass density) 


2. Constraints: 


u z = 6 x = 0 y = 0 (all points) 
u y = 0 (point 101 ) 
u x = u y = 0 (point 1) 

3. Loads. 


F icn ,x " 3,423,17 

B =1.0 (default load factor) 


C. Results 


The theoretical natural frequency for the first mode is given by 


f = 


1 


4pA Z 2 


(^-f) 


- Va 


Hertz 


For this loading of one half the Euler buckling value, the theoretical value is 14.6269 Hertz 
for the bending mode. The NASTRAN result is 14.62325 Hertz. 
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Figure 1- One hundred cell beam. 
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RIGID FORMAT No. 14, Static Analysis with Cyclic Symmetry 
Circular Plate Using Cyclic Symmetry (14-1-1) 


A. Description 

A constant thickness circular plate with six radial stiffeners and a central hole, shown in 
Figure 1, is analyzed using dihedral symmetry. The plate is subjected to a uniform pressure load 
applied over a 60° segment of the plate. 

The finite element model is shown in Fiqure 2. The stringers are 60° apart but only 30° of the 
structure needs to be modeled when using the dihedral symmetry option. There are 12 subcases since 
these are 2 half segments in a 60° segment and only one loading condition. The CYJ0IN bulk data 
card defines those points in the middle of the segment (SIDE 2) and those points on the boundary 
between segments (SIDE 1). 

B. Input 

1. Parameters. 


R = 1.0 

0 

(outside radius) 

R = .14 

i 

(inside radius) 

t = .01 

(plate thickness) 

a = .06 

(height and width of stiffeners) 

E = 10. 6x1 U 6 

(modulus of elasticity) 

v = .325 

(Poisson's ratio) 


2. Boundary Conditions: 

U r = U Q = 6 z = 0 (all points) 

U z = 6 r = 0 (along r = 1.0) 

3. Appl ied loads : 

Pressure = 200.0 between 0 = 60° and 120° 

4. Cyclic symmetry parameters: 

CTYPE = DRL 

KMAX = 2 
NSEGS = 6 
N10AD = 1 
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C. Results 


The structure can be analyzed using rotational symmetry or dihedral symmetry described here 
and the results will be identical. 

The results for the normal displacements are given in Table 1 for r = 0.46. 


/ 


14.1-2 (12/31/74) 






Pressure load area 



i 


T 

.01 


=J I j -06 

I I -025 ♦ 

Section A-A 


Figure 1. Circular plate with stiffeners 
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0.46 


SIDE 1 


Figure 2 



Finite element model. 
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RIGID FORMAT No. 15, Normal Modes Analysis Using Cyclic Symmetry 
Modal Analysis of a Circular Plate Using Cyclic Symmetry (15-1-1) 


A. Description 

The natural frequencies of a constant thickness circular plate with six radial stiffeners and 
a central hole are obtained using the rotational symmetry option. The structure, shown in Figure 
1, is simply supported at the outer circumference. 

The finite element model is shown in Figure 2 representing only sixty degrees of the plate. 
Note that since the stiffeners are on the symnetry boundary, only 1/2 of the actual properties 
are used. The bulk data cards demonstrated are the CYJ0IN and PARAM. 


B. Input 


1. Parameters: 


R o 

R i 

t 

a 

E 

v 

P 


1.0 

(outside radius) 

.14 

(inside radius) 

.01 

(plate thickness) 

.06 

(height and width of stiffeners) 

10. 6x7 0 6 

(modulus of elasticity) 

.325 

(Poisson's ratio) 

2.59xl0- 4 

(mass density of plate and stiffeners) 


2. Boundary conditions: 


u r = u Q = 6 z = 0 (all points) 
u z = 0 r = 0 (along r = 1 .0) 

3. Eigenvalue extraction data. 


Method: Inverse power 

Region of interest: 0.0 < f < 8000 
Number of desired roots: 3 
Normalization: maximum 


4. Cyclic symmetry parameters: 

CTYPE R0T 
KINDEX 2 
NSEGS 6 
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C. Results 


Solutions can be obtained using the dihedral symmetry or rotational symmetry described here. 
Results are accurate to approximately six significant figures. 


Table 1. Natural Frequencies 


Mode 

Frequency 

i 

4288.2 

2 

6844.3 

3 

11524.3 
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Figure 2. Finite element model. 
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